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ABSTRACT 


This  Second  Quarterly  Technical  Report  for  the  LASA  Experimental 
Signal  Processing  System  discusses  the  effort  expended  during  the  second 
quarter  to  provide  the  hardware  and  software  to  support  research  and 
development  directed  toward  the  study  of  seismic  signal  processing;  and 
it  delineates  work  tasks  planned  for  the  next  quarter.  This  document  also 
presents  detailed  information  related  to  configurations,  array  monitor  and 
control,  data  compression,  scaling,  display,  and  array  design. 
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INTRODUCTION 

The  work  reported  herein  was  performed  under  the  "LASA  Experimental 
Signal  Processing  System,"  Contract  Number  F19628-67-C-0198,  and  is  a 
continuation  of  the  "Large  Aperture  Seismic  Array  Signal  Processing  Study,"  1 
Contract  Number  SD-296,  and  the  "Large  Aperture  Seismic  Array  Signal  Pro¬ 
cessing  Communications  and  Simulation  Study,"2>3>4  Contract  Number 
AF19(628)-5948. 

The  purpose  of  this  contract  is  to  design,  develop  and  implement  the  LASA 
Experimental  Signal  Processing  System  (ESPS),  including  the  hardware  and 
software,  to  provide  an  experimental  capability  to 

a.  Evaluate  performance  of  the  system  in  accordance  with  system 
requirements. 

b.  Demonstrate  the  capability  to  meet  basic  LASA  signal  processing 
objectives . 

c.  Conduct  research  to  develop  means  of  improving  and  extending  the 
capability  of  the  system. 

d.  Perform  seismic  and  signal  processing  experiments  of  interest. 

In  addition  to  the  above,  activities  to  support  the  future  integration  of  the 
system  will: 

a.  Specify  the  requirements  to  interface  with,  control  and  monitor  the 
LASA  array. 

b.  Specify  the  requirements  for  a  facility  to  house  the  system  in  Montana. 
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During  this  second  reporting  period,  effort  was  principally  directed  to 
three  areas: 

1.  Completion  of  plans  and  occupancy  of  the  Seismic  Array 
Analysis  Center  (SAAC),  Van  Ness  Centre,  Suite  450,  4301 
Connecticut  Avenue,  Washington,  D.  C.  20008. 

2.  Installation  and  utilization  of  the  Detection  and  Event  Processors 
and  peripheral  equipment  as  System/360  Model  40  G  and  H 
general-purpose  computers,  respectively. 

3.  Initiation  of  array-processing  synthesis,  analysis  and  design 
for  Array  II. 


1-2 


Section  2 


RESULTS  AND  SUMMARY  OF  WORK  ACCOMPLISHED 

The  three  subsections  which  follow  correspond  to  the  functional  disciplines. 

2.1  SYSTEM  ACTIVITIES 

Studies  during  this  quarter  addressed  the  composite  detection  system  and 
developed  a  comprehensive  array-processing  synthesis  procedure  responsive 
to  both  seismic  missions  of  event  detection  and  location.  Analysis  of  array 
geometry,  communications,  process  errors,  and  implicit  computer  burden  con¬ 
sidered  the  inherent  data  character;  these  analyses  are  presented  in  the 
appendixes. 

Since  large  seismic  arrays  produce  extensive  raw  data,  message  encoding 
is  a  critical  system  constituent.  Delta  modulation,  sample  rate  reduction,  and 
floating-point  word  formats  have  been  examined  to  establish  data  compression 
techniques  appropriate  to  array  data  acquisition  and  communication  (see  Appen¬ 
dix  C).  We  recommend  use  of  a  base  4  modified  floating-point  encoding  of  short 
period  seismic  data  sampled  at  10  Hz.  This  study  considered  system  perform¬ 
ance  implications,  subarray  equipment  compatibility,  and  the  objective  to  effi¬ 
ciently  use  commercially  available  modular  communication  facilities. 

The  resulting  communication  rate  is  thereby  reduced  to  approximately 
35  percent  of  current  requirements  without  significantly  impairing  system  per¬ 
formance.  Actual  seismic  histories  were  employed  to  evaluate  several  trace 
reconstruction  algorithms,  and  the  observed  errors  were  commensurate  with 
classical  error-bound  estimates.  The  introduced  noise  power  and  deviation 
bounds  are  negligible  when  reconstruction  is  derived  from  four  or  more  data 
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points,  which  verifies  the  signal  integrity  at  this  reduced  sampling  rate.  Instru¬ 
ment  rescaling  and  analog  filter  tuning  adjustments  are  latent  modifications 
which  can  further  improve  system  performance. 

Detection  process  parameters  have  been  reviewed,  and  an  initial  scaling 
has  been  formulated  which  is  compatible  with  operational  developments,  evalu¬ 
ation,  and  data  analysis  support.  Sampling  rate,  data  quantum,  saturation  level, 
and  estimated  statistical  characteristics  of  the  composite  noise  inherent  in  both 
the  raw  data  and  the  fixed-point  digital  process  have  been  evaluated  at  the  nodes 
of  each  stage  (see  Appendix  D).  Numerical  error  contributions  are  insignificant 
thereby  substantiating  the  detection  architecture.  Filter  response  deviations 
due  to  coefficient  quantization  are  presented,  and  a  slight  pass  band  ripple  is 
evident  for  conventionally  rounded  values.  The  encountered  characteristics  are 
adequate,  and  since  the  detection  filter  can  be  subsequently  upgraded  when 
statistically  significant  spectral  estimate  ensembles  are  available,  no  secondary 
filter  adjustments  were  considered  necessary  at  this  time. 

Experimental  display  modification  to  support  a  strip  chart  recorder  proved 
feasible  and  is  described  in  Appendix  E.  This  extension  implements  a  rapid 
multichannel  hard-copy  background  capability  of  medium  accuracy  which  facili¬ 
tates  process  documentation  and  provides  an  efficient  method  for  gross  data 
editing. 

Transformations  between  the  frequency  and  time  domains  are  often  desired 
to  efficiently  analyze  signals  and  characterize  processes.  Array  instrument 
calibration,  phase  velocity  estimation,  reverberation  suppression,  and  special 
event  processing  are  operational  functions  which  may  conceivably  require  the 
application  of  the  Cooley-Tukey  (Fast  Fourier)  algorithm.  An  analysis  of  the 
errors  encountered  when  implementing  this  process  with  16-bit  fixed-point 
arithmetic  was  conducted  and  simulations  performed.  When  proper  inter¬ 
mediate  scaling  precautions  are  exercised,  the  concept  appears  feasible. 
Estimates  of  implementing  this  operation  as  a  microcoded  algorithm  have 
now  been  generated,  and  an  execution  time  of  100  milliseconds  as  opposed 
to  920  milliseconds  in  fixed-point  and  1080  milliseconds  in  floating-point,  coded 
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in  assembly  language,  are  required  to  transform  a  256  complex  number  array 
on  an  IBM  360/2040.  As  the  Event  Processor  architecture  evolves,  significance 
of  this  advancement  to  the  system  will  be  evaluated  to  determine  the  desirability 
of  implementing  the  transformation  as  an  addition  to  the  microcode  repertoire. 

The  comprehensive  array  synthesis  presented  in  Appendix  F  illustrates  a 
methodical  procedure  which  is  pertinent  to  overall  system  objectives,  provides 
tradeoff  flexibility,  and  identifies  crucial  variables  which  require  measurements 
to  establish  experimental  verification.  This  establishes  the  necessary  objectives 
of  a  precursor  data  measurement  program  and  provides  the  structure  for  sub¬ 
sequent  array  designs. 

2.2  PROGRAMMING 

Programming  effort  during  this  quarter  centered  on 

1.  Microprogramming 

2.  Detection  Processing 

3.  Supervisor/Monitor 

4.  Program  Conversion. 

Accomplishments  in  each  of  these  areas  is  presented  in  the  following  subsections. 
2.2.1  Microprogramming 

Microprogramming  activity  centered  around  two  areas:  continued  testing 
and  the  new  RIT  threshold  computation. 

Following  installation  of  the  2040G  at  SAAC,  microcode  testing  was  con¬ 
tinued.  After  completion  of  these  tests,  debug  work  on  the  detection  controller 
revealed  a  number  of  data-dependent  errors;  these  were  corrected  as  they 
appeared  with  no  loss  in  execution  speed. 

The  improved  RIT  threshold  computation  operation  (previously  reported 
upon5)  was  installed  in  both  machines.  Testing  of  this  addition  revealed  no 
operational  errors. 
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2.2.2  Detection  Processing 

Detection  processing  activity  was  centered  in  two  areas:  data  acquisition 
and  the  detection  controller. 

2. 2.2.1  Data  Acquisition 

The  initial  version  of  this  program5  was  modified  to  take  advantage  of  the 
real-time  input  format.  The  program,  operating  in  an  off-line  mode,  accepts 
as  input  LASA  data  tapes  and  generates  a  seismometer  level  tape  in  the  same 
format  as  the  real-time  seismometer  input  format.  This  data  tape  will  serve 
as  the  primary  input  to  the  Detection  Processor.  The  editing  and  reformatting 
function  will  also  provide  a  record  of  the  type  and  number  of  errors  (modem, 
parity  and  illegal  values)  on  an  individual  seismometer,  subarray  and  array 
basis. 

2. 2. 2. 2  Detection  Controller 

Work  was  aimed  at  expanding  the  version  written  for  microcode  testing. 
Initial  operating  parameters  were  defined  as  were  input  and  output  formats.  A 
version  of  the  controller  was  completed  and  initial  experimentation  was  begun. 

The  operating  parameters  are  previously  reported.4  These  parameters 
were  used  to  assess  the  operating  posture  of  the  controller.  Input  data  consisted 
primarily  of  a  special  data  tape  which  contained  six  known  events.  A  pseudo¬ 
noise  tape  was  also  created  and  used  to  obtain  an  initial  noise  estimate  for 
insertion  into  the  noise  averaging  section.  The  type  of  outputs  generated  are 
also  previously  reported.4 

Testing  of  the  controller  pointed  up  some  errors  in  the  microcodes,  and 
these  in  turn  resulted  in  some  reprogramming  effort.  In  particular,  core 
storage  allocation  and  arrangement  were  the  major  areas  of  redefinition. 

Initial  experimentation  consisted  of  running  the  detection  controller  to 
generate  intermediate  data  tapes  of  subarray  beams  and  LASA  beams.  These 
tapes  were  plotted  and  the  plots  used  to  locate  errors  in  the  controller.  The 
listed  output  of  detected  arrivals  was  also  an  aid  in  error  detection. 
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2.2.3  Supervisor/Monitor 

Upon  arrival  of  the  IBM  System/360  Model  2040G  at  SAAC,  significant  effort 
was  expended  in  systems  generation  and  maintenance  functions.  The  operating 
system  in  use  is  a  modified  version  of  the  IBM  Disk  Operating  System  (DOS) 
under  which  the  ESPS  programming  development  will  proceed.  Supervisor/ 
monitor  activity  for  the  Detection  and  Event  Processors  is  discussed  in  the 
following  sections. 

2. 2. 3.1  Detection  Subsystem 

Performance  Test  Tool.  A  means  of  measuring  the  detection  processing 
capability  under  conditions  of  minimum  overhead  was  designed  and  coded.  The 
object  was  to  operate  independently  of  DOS  so  as  to  maximize  the  available 
computing  time.  This  allowed  a  measure  of  the  upper  bound  of  computing  capa¬ 
bility.  Also,  by  comparison  with  a  detection  processing  program  which  used 
DOS,  a  measure  of  the  overhead  incurred  in  the  latter  arrangement  would  be 
available.  This  would  aid  in  the  planning  of  the  eventual  real-time  detection 
monitor.  This  test  tool  was  made  available  so  it  could  be  combined  with  the 
detection  processing  program. 

Prototype  Detection  Monitor.  Work  progressed  in  the  development  of  the 
detection  monitor.  The  first  version  will  be  capable  of  taking  simulated  real¬ 
time  inputs  asynchronously  with  detection  processing.  This  requires  a  multi¬ 
programming  capability.  The  modifications  to  the  DOS  needed  to  make  this 
possible  were  identified. 

2. 2. 3. 2  Event  Subsystem 

The  Disk  Operating  System  does  not  fully  support  the  implementation  of 
programs  using  the  1627  plotter,  the  experimental  display,  channel-to-channel, 
and  the  2250  display.  The  provision  of  a  general-purpose  capability  tailored  to 
the  SAAC  configuration  was  given  the  most  emphasis  during  this  quarter.  Speci¬ 
fications  were  made  for  the  needed  supervisor  modifications  to  support  the 
plotter  and  the  experimental  display  on  the  2040H. 
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The  first  objective  in  developing  the  experimental  Event  Processor  monitor 
includes  as  a  fundamental  requirement  the  ability  to  execute  programs  auto¬ 
matically.  Whereas  in  a  general-purpose  operating  system,  the  load  on  the 
system  is  controlled  by  a  computer  operator,  in  an  automatic  system  a  job  is 
executed  by  signals  set  and  detected  internally.  For  example,  event  processing 
will  be  triggered  by  an  interrupt  via  the  channel -to-channel  from  the  Detection 
Processor  subsystem.  The  modifications  to  DOS  needed  to  build  this  capability 
were  investigated  during  this  quarter. 

The  capability  represented  by  the  prototype  system  would  be  extended  in 
further  versions  to  provide  the  automatic  job  scheduling  capability  to  build  a 
system  responsive  to  its  external  environment. 

2.2.4  Program  Conversion 

The  programs  developed  for  the  IBM  7090  during  previous  studies  were 
converted  to  the  IBM  System/360  to  run  under  the  IBM  Disk  Operating  System. 
In  general,  this  conversion  process  was  easily  accomplished  since  the  majority 
of  the  programs  were  written  in  FORTRAN.  The  major  changes  made  to  the 
FORTRAN  programs  involved  array  sizes  and  logical  statements. 

Two  programs,  the  LASA  tape  edit  and  the  read  routine,  were  rewritten. 
The  LASA  tape  edit  program,  written  in  Assembly  Language  for  the  IBM  7090, 
finds  its  System/360  counterpart  in  the  data  acquisition  program  discussed 
above.  The  read  routine  provides  the  capability  for  the  FORTRAN-oriented 
data  analysis  programs  to  process  the  data  tapes  prepared  by  the  data  acqui¬ 
sition  program.  It  is  an  input  routine  which  accepts  the  real-time  formatted 
tapes,  and  unpacks  and  places  in  floating-point  form  the  seismometer  channels 
on  these  tapes. 

The  conversion  process  offered  the  opportunity  to  increase  the  capability 
of  several  programs.  In  particular,  program  capability  in  the  signal  and  noise 
power  calculations  and  correlation  areas  was  extended.  The  power  programs 
were  provided  with  a  sliding  integration  window  capability  as  opposed  to  the 
fixed  window  of  the  IBM  7090  version.  In  addition,  power  averages  for  25,  19, 
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16  and  13  seismometers,  and  21,  17,  13  and  9  subarrays  were  included  in  the 
calculations. 

The  correlation  program  was  modified  to  correlate  pairs  of  data  channels, 
provide  intermediate  correlation  coefficients  for  all  pairs,  and  record  the  output 
information  as  a  function  of  distance  between  the  sensors  in  each  pair.  The 
output  format  includes  a  running  sum  of  the  correlation  coefficients  for  all  pairs 
considered  in  the  processing  from  which  an  average  correlation  coefficient  as 
a  function  of  distance  between  elements  can  be  obtained. 

Minor  modifications  were  made  to  several  other  programs  to  provide 
greater  flexibility  and  efficiency  in  their  execution  capability. 

2.3  OPERATIONS 

Effort  during  this  quarter  was  expended  in  the  following  operational  areas: 

1.  SAAC 

2.  Maintenance  Plan 

3.  System  Procedures 

4.  Experimental  Display 

5.  ESPS  General  Specification 

6.  ESPS  Acceptance  Test  Plan 

7.  Interface  Monitor  and  Control  Specification. 

The  activity  in  each  of  these  areas  is  discussed  in  the  following  sections. 
2.3.1  SAAC 

The  Seismic  Array  Analysis  Center  (SAAC)  was  occupied  by  IBM  LASA 
Project  personnel  on  March  20,  1967.  Prior  to  this,  specific  office  layouts  and 
the  computer  room  layout  were  prepared  and  approved;  and  coordination  of  the 
computer  installation  commenced. 

Installation  of  the  initial  SAAC  machine  system  (described  in  Appendix  A) 
was  completed,  and  the  GSA  acceptance  requirements  were  met  during  March 
1967. 
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As  program  debug  efforts  began,  critical  examination  was  given  to  system 
utilization  and  performance.  As  a  result  and  in  view  of  current  overall  project 
needs,  several  changes  were  made  to  adapt  the  hardware  system  to  a  more 
efficient  configuration  for  software  development.  These  changes  include: 

a.  System  interconnections  will  be  implemented  on  an  as-required 
basis  to  minimize  software  interaction. 

b.  Additional  unit-record  equipment  (line  printer,  reader/punch)  was 
ordered  to  provide  general-purpose  capabilities  on  the  System/360 
2040G  (Detection  Processor)  during  the  programming  debug  phase. 

c.  The  conflicts  arising  from  use  of  a  shared  disk  control  unit,  when 
the  two  processors  are  operating  independently,  results  in  inefficient 
processor  utilization.  Since  this  will  be  a  common  operational  mode 
during  system  development,  an  additional  disk  control  unit  and  drive 
were  ordered. 

d.  Plans  were  formed  to  establish  the  System/360  Model  2044H  as  a 
stand-alone  computing  system  with  the  capability  to  be  integrated 
as  an  interconnected  processor  at  a  later  date. 

The  Detection  Processor  was  used  for  approximately  200  hours  for  this 
quarter,  with  the  major  amount  of  this  time  being  used  for  debugging  programs 
and  for  data  analysis. 


2.3.2  Maintenance  Plan 

During  the  reporting  period,  attention  was  given  to  developing  a  maintenance 
approach  for  the  system  hardware  tailored  to  the  configuration  being  used  and 
the  overall  system  performance  requirements. 

Basic  maintenance  concepts  are  valid  in  the  usual  situation  where  the 
processor  can  be  reconfigured  to  maintenance  subsystems  allowing  the  prime 
function  of  detection  processing  to  proceed  uninterrupted.  An  unusual  case 
results  when  "machine  state  failures"  (those  which  occur  when  a  precise  pro¬ 
gramming  system  environment  has  been  established)  occur  between  the  inter¬ 
connected  processors  and  their  associated  shared  control  units.  Several  tech¬ 
niques  have  evolved  in  other  systems  of  the  same  general  type  and  these  are 
being  studied  to  determine  the  best  approach. 
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2.3.3  System  Procedures 


Two  undesirable  situations  were  identified  as  a  result  of  investigating  the 
fail  soft  technique  to  be  employed  in  the  on-line,  real-time  system.  These  were: 

a.  Dominant  control  unit  —  the  last  control  unit  of  a  series  of  control 
units  provides  the  termination  for  the  channel/control  unit  signals. 

When  power  to  this  unit  is  removed,  or  is  lost,  a  floating  condition 
occurs  on  these  lines  which  can  render  a  processor  inoperative. 

Due  to  the  interconnected  configuration  in  the  system,  this  condition 
is  intolerable.  The  designed  solution  will  eliminate  the  requirement 
for  power-up  on  the  terminating  control  unit. 

b.  Power  on/ off  transients  —  during  the  power  on/off  sequences  of 
control  units,  undesirable  noise  signals  occur  on  the  interface  con¬ 
trol  and  signal  lines.  As  a  result,  error  conditions  can  be  generated 
which  would  adversely  affect  system  operation.  As  in  item  a.,  the 
presence  of  the  system  interconnections  makes  this  an  unacceptable 
situation.  A  design  to  isolate  any  unit  w'hich  is  transitioning  power 
from  interface  control  and  signal  lines  until  transition  has  been 
completed  is  in  progress.  The  required  equipment  changes  have 
been  requested  for  installation  as  soon  as  the  field  change  packages 
are  available. 


2.3.4  Experimental  Display 

The  logic  design  and  fabrication  of  the  hardware  for  the  on-line  experimental 
display  system  has  been  completed.  The  system  is  now  under  unit  test.  Design 
reviews  have  been  conducted  and  the  necessary  documentation  for  software 
generation  has  been  provided. 


2.3.5  ESPS  General  Specification 

A  preliminary  version  of  the  ESPS  General  Specification  was  written  and 
distributed  for  project  review.  This  specification  establishes  performance, 
design,  development  and  test  requirements  for  the  LASA  Experimental  Signal 
Processing  System.  Equipment  selection  and  program  design  necessary  to  pro¬ 
vide  a  system  capable  of  supporting  research  and  development  in  techniques  of 
seismic  signal  processing  are  presented. 
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This  document  provides  the  inputs  which  form  the  basis  for  developing 
further  system  documentation  such  as  Programming  Specifications  and  the 
Acceptance  Test  Plan.  In  turn,  these  documents  will  provide  the  criteria  for 
design,  performance  and  test  of  the  system,  and  system  components. 

2.3.6  ESPS  Acceptance  lest  Plan 

The  initial  draft  of  the  Acceptance  Test  Plan  was  completed  and  distributed 
for  project  review  during  this  reporting  period.  This  document  presents  the 
tests,  equipment,  program,  and  subsystems  necessary  to  demonstrate  the  capa¬ 
bility  of  the  system  to  meet  the  specified  requirements. 

Included  in  this  document  are  system  descriptions,  test  requirements,  con¬ 
cept  of  testing,  test  methods  and  test  support  requirements.  Information  supplied 
by  this  plan  will  be  used  as  a  guide  in  developing  the  quality  assurance  portion 
of  the  programming  specifications.  This  portion  of  the  programming  specifica¬ 
tion  will  establish  the  details  governing  each  of  the  tests  to  be  conducted. 

2.3.7  Interface  Monitor  and  Control 

The  interface  monitor  and  control  requirements  (Appendix  B)  were  com¬ 
pleted  during  this  reporting  period.  The  preliminary  version  of  the  Array  Inter¬ 
face,  Monitor,  and  Control  Functional  Specification  was  drafted  and  has  been  sub¬ 
mitted  for  project  review  prior  to  final  editing  and  release. 
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Section  3 


PLANS 

The  following  activities  are  planned  for  next  quarter.  Discussions  correspond 
to  the  functional  disciplines.  Each  task  is  presented  with  sufficient  perspective  to 
identify  its  interrelation  and  contributions  to  the  project. 

3.1  SYSTEM  ACTIVITIES 

Threshold  parameters  will  be  studied  to  characterize  the  noise  rejection  and 
signal  identification  capabilities  of  the  implemented  detection  criterion.  An 
initial  set  of  parameters  will  be  selected  to  complete  the  microcode  definition  for 
data  analysis. 

A  dynamic  array  geometry  model  will  be  developed  to  establish  the  optimum 
gain  array  configuration  within  a  circular  constraint  for  specific  homogeneous 
noise  fields.  This  model  will  support  a  performance/cost  tradeoff,  and  can 
potentially  extend  the  art  of  prudent  array  design.  Noise  field  models  are  an 
essential  input  to  this  activity,  and  the  capability  to  derive  such  empirical 
characterizations  will  be  actively  pursued. 

Comprehensive  data  analysis  activities  will  be  structured  to  efficiently  provide 
additional  base  data  for  phase  data  anomaly  characterization  and  to  extend  the 
knowledge  of  both  the  signal  parameters  and  the  inherent  noise  field.  Subsequent 
process  evaluation  objectives  will  be  considered  to  minimize  the  composite  data 
and  computer  requirements. 

Automatic  phase  delay  estimation  techniques  will  be  pursued  in  an  effort  to 
develop  a  precise  analytical  capability  amenable  to  the  efficient  acquisition  of 
statistically  significant  anomalous  observations  from  the  available  bulk  raw  data. 
Samples  obtained  from  the  presently  defined  regions  or  sectors  will  be  used  to 
evaluate  the  anomaly  characterization  prior  to  being  incorporated  into  the  library 
data. 
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Phase  delay  retrieval  and  signal  wavefront  analysis  segments  of  the  event 
process  structure  will  be  detailed  to  permit  the  initiation  of  precursor  encoding 
for  technique  evaluation.  Initial  linkage  requirements  will  thus  be  identified,  and 
will  serve  as  departure  datums  for  the  evaluation  of  the  event  process. 

A  standardized  orthogonal  polynomial  representation  will  be  adopted  for  the 
teleseismic  range.  Error  characteristics  will  be  examined  for  several  event 
arrival  phases  to  determine  if  single  precision  computation  is  adequate.  Geographic 
coordinate  estimation  will  be  examined  to  identify  the  necessary  program  structure. 

Beam  pattern  shaping  will  be  examined  to  determine  if  such  techniques  are 
effective  in  combating  coherent  noise  originating  in  the  proximity  of  desired 
signals.  If  feasible,  this  would  significantly  extend  the  special  event  processing 
capability  to  enhance  specific  signals  in  the  presence  of  simultaneous  wave 
arrivals. 

3.2  PROGRAMMING 

Plans  for  the  next  quarter  include  activities  in  both  the  Detection  and  Event 
Processor  functions. 

3.2.1  Detection  Processor 

The  detection  controller  routine  will  be  modified  so  that  it  will  accept  as  input 
the  edited  and  reformatted  tapes  produced  by  the  data  acquisition  program.  An 
additional  version  of  this  routine  will  be  produced  and  will  contain  the  required 
modifications  to  allow  its  use  in  data  analysis  activity. 

The  first  model  of  the  detection  monitor  will  be  generated  and  tested  by 
simulation  of  the  detection  controller  and  simulation  of  a  real-time  interface  using 
the  event  processor  and  channel -to-channel  capability.  Following  the  successful 
testing  of  this  model  it  will  be  combined  with  the  detection  controller  and  will  form 
the  first  model  of  the  Detection  Processor  program. 
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3.2.2  Event  Processor 


The  required  modifications  and  additions  to  the  Disk  Operating  System  will 
be  generated  and  will  provide  the  required  support  for  the  special  peripheral 
input/output  devices  attached  to  theEvent  Processor  (2040H).  This  support  includes 
the  1627  plotter,  experimental  display,  channel -to-channel,  and  the  2250  visual 
display  unit.  This  activity  will  provide  a  general-purpose  operating  system  tailored 
to  the  SAAC  configuration. 

Generation  of  a  multiprogramming  capability  will  provide  efficient  and 
flexible  use  of  the  machine  configuration.  The  presence  of  two  computers  with 
one  printer,  one  plotter  and  one  card  punch  necessitates  the  spooling  of  output 
from  tape  to  these  devices. 

The  IBM  Disk  Operating  System  with  multiprogramming  permits  the  con¬ 
current  operation  of  three  independent  jobs.  To  facilitate  the  spooling  of  output, 
we  will  provide  utility  programs  which  will  operate  full-time  in  two  of  these 
partitions,  while  the  third  partition  will  have  full-stacked  job  capability. 

Experimental  program  development  will  be  effected  to  aid  in  developing  Event 
Processor  functions.  This  program  will  include  activity  in  the  time  delay,  beam 
selection,  arrival  identification,  and  event  location  areas. 

To  facilitiate  the  use  of  the  1627  plotter,  a  set  of  basic  plotter  routines  will 
be  developed.  A  partial  plotter  routine  program  package  is  available,  and  the 
required  modifications  and  extensions  to  this  package  will  be  made  to  provide  the 
plotter  capability  required  at  SAAC. 

The  programs  required  to  present  waveforms  and  beams  on  the  experimental 
display  will  be  developed  and  tested  when  this  display  is  integrated  with  the  2040H 
processor.  The  definition  and  design  of  the  program  specifications  for  the  2250 
visual  display  unit  will  be  addressed,  and  implementation  of  these  functions  will 
be  in  progress. 

3.3  OPERATIONS 

Plans  for  the  next  quarter  include  acitivity  in  SAAC  expansion,  equipment, 
and  data  item  reporting. 
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Tentative  office  layout  plans  have  been  drafted  to  expand  SAAC  approximately 
2800  square  feet.  These  plans  include  the  addition  of  a  seminar  room,  several 
offices,  a  small  library  and  document  control  area.  Other  modifications  of  office 
space  will  result  in  changing  the  location  of  the  reception  area  and  the  possible 
addition  of  one  or  two  more  offices  by  conversion  of  one  conference  room.  Plans 
for  rearrangement  of  the  computer  rooms  and  tape  storage  area  are  in  progress. 
These  plans  will  consider  alternative  configurations  of  the  interface  system 
equipment.  Operation  of  SAAC  has  demonstrated  the  need  for  more  storage  space 
for  data  processing  supplies  and  magnetic  tapes  as  well  as  desk  space  for  operation 
of  the  tape  library. 

The  installation  of  the  LASA  intermediate  computing  system  is  scheduled 
for  completion  at  the  start  of  the  fourth  quarter.  Heavy  activity  is  expected  in 
integrating  the  supervisor  software  routines  required  to  support  the  various  system 
interconnections.  With  both  the  System /360  40G  and  System /360  40H  fully 
operational  during  the  next  quarter,  an  expected  partial  second  shift  will  be  operating 
to  fully  use  the  176  hours  of  machine  time  allowed  in  the  basic  monthly  rental.  A 
large  shipment  of  tapes  will  be  ordered  for  processing  through  our  edit  programs. 

In  the  future,  tape  ordering  will  be  on  a  continuing  basis  in  order  to  develop  an 
adequate  base  for  system  experimentation  and  test. 

The  LASA  experimental  display  system  modified  for  on-line  operation  is 
scheduled  for  hardware  integration  with  the  LASA  computing  system  during  the 
third  quarter.  Primary  emphasis  will  be  placed  on  the  beam  and  waveform 
functions  so  that  the  program  debug  effort  can  be  started.  The  strip  chart  recorder 
functions  will  be  integrated  when  test  software  is  available. 

Immediate  plans  call  for  project  review  of  the  General  Specification  and 
Acceptance  Test  Plan  so  that  the  next  level  of  system  documentation  may  in  turn 
be  initiated. 
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Appendix  A 

MACHINE  CONFIGURATIONS 

The  initial  and  second  SAAC  machine  configurations  are  presented  in 
figures  A-l  and  A-2  respectively. 

During  March  1967  the  IBM  System/360  Model  40G  (Detection  Processor) 
was  received  and  installed,  and  on  March  20,  1967,  it  was  released  to  the  LASA 
Project  for  utilization.  In  April  1967,  the  IBM  System/360  Model  40H  (Event 
Processor)  was  received  and  installed.  It  is  expected  that  it  will  be  accepted 
for  use  in  May  1967. 

Figure  A-2  presents  the  SAAC  machine  configuration  as  it  will  exist  during 
the  next  quarter.  The  cross-hatching  indicates  units  to  be  added  to  the  system 
during  the  quarter. 


A-l 


A-2 


Figure  A-l.  Initial  Configuration 
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Figure  A-2.  Second  Configuration 


Appendix  B 

ARRAY  MONITOR  AND  CONTROL  REQUIREMENTS 
B.l  INTRODUCTION 

An  operator's  console  and  a  maintenance  console  shall  be  provided  to  facili¬ 
tate  monitoring  the  seismometer  outputs,  as  well  as  the  computer  processing, 
to  make  printed  records  of  the  signals  being  processed,  and  to  initiate  commands 
to  the  processor  or  array  to  perform  diagnostic  tests  for  maintenance  purposes 
or  to  calibrate  the  array  seismometers.  Operation  of  the  console  displays  and 
keyboard  entries  shall  be  in  conjunction  with  computer  programs  which  control 
the  transfer  of  data  between  the  computer  and  the  console. 

B.2  GENERAL  DESCRIPTION 

B.2.1  Array  Monitoring 

The  rectified  output  of  each  seismometer  shall  be  integrated  over  a  specified 
time  interval.  Selected  seismometer  outputs  may  be  monitored  by  displaying 
the  signal  received  on  a  cathode  ray  tube  waveform  display.  Six  to  eight  such 
signals  may  be  displayed  simultaneously  for  comparison,  and  any  particular 
signal  shall  be  selectable  to  be  displayed  continuously  for  detailed  examination. 

In  addition  to  the  waveform  display,  a  capability  must  exist  for  making  a  printed 
record  of  the  signals  on  a  strip  chart  recorder  or  similar  analog  recording 
device  for  later  study  or  for  examination  at  a  point  remote  from  the  processing 
center. 
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Array  alarm  conditions  shall  be  continuously  monitored  to  provide  indication 
of  power  failure,  etc.  Since  the  operation  of  the  array  will  be  dependent  upon 
the  performance  of  the  two  way  communication  link  between  the  array  and  the 
processing  center,  it  will  be  necessary  to  monitor  the  communication  equipment 
and  to  provide  indicators  to  the  operator,  as  well  as  to  the  processor,  in  the 
event  of  a  communication  link  failure.  The  error  detection  bits  of  data  received, 
modem  AGC  warning,  and  loss  of  synchronism  of  message  frame  must  be 
monitored. 

B.2.2  Monitor  Computer  Processing 

Computer  processing  will  be  monitored  on  a  waveform  display  and/or  a 
beam  display.  These  displays  will  make  it  possible  to  view  a  large  quantity 
of  output  data  as  a  single  image  and  to  obtain  information  on  absolute  as  well 
as  relative  values. 

The  computer  equipment  shall  be  monitored  and  status  indicated. 

B.2.3  Array  Control  Commands 

A  number  of  test  facilities  have  been  designed  into  the  subarray  equipment 
and  will  be  placed  in  operation  by  commands  from  the  computer  or  consoles. 
These  facilities  are  necessary  for  maintenance  of  the  array,  as  well  as  calibra¬ 
tion  of  the  array  seismometers,  and  provide  a  means  of  sampling  special  sensing 
equipment  at  the  subarrays.  A  designated  six -bit  field  of  the  message  trans¬ 
mitted  to  each  subarray  shall  contain  the  command  codes.  These  command 
codes  will  be  transmitted  to  any  one  selected  subarray  at  a  time  and  must  not 
interfere  with  normal  operation  of  the  other  subarrays.  The  formation  of  the 
commands  and  the  selection  of  the  subarray  and  subarray  data  must  be  per¬ 
formed  under  program  control  or  by  manual  entry  keys  on  the  consoles.  The 
command  code  received  at  the  subarray  shall  be  transmitted  back  to  the  Data 
Center  for  verification. 
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B.2.4  Output  Display  for  Maintenance  of  Array  and  Communication  Link 

One  waveform  display  will  be  used  for  monitoring  the  array  and  processor 
and  a  second  waveform  will  be  used  for  maintenance  purposes.  A  strip  chart 
recorder  will  provide  a  permanent  record  of  the  output  from  any  seismometer. 
Most  of  the  diagnostic  routines  will  be  performed  automatically,  under  program 
control,  with  their  results  summarized  and  a  printed  console  log  presented  on 
the  output  printer.  When  desirable,  the  spectral  response  of  a  selected  channel 
may  be  displayed  on  the  waveform  display  with  acceptable  limits  superimposed. 
Or,  the  numeric  evaluation  of  the  spectral  response,  as  well  as  numeric  parame¬ 
ters  or  status  indicator,  may  be  displayed. 

B.3  OPERATIONAL  CHARACTERISTICS 

B.3.1  Digital  to  Analog  Conversion 

Eight  digital  to  analog  conversion  units  shall  be  provided  to  give  eight 
independent  analog  output  channels.  These  channels  will  accept  digital  informa¬ 
tion  in  binary  form  (13  bits  plus  sign  and  provide  analog  output  voltage  of  +  10 
volts  with  an  accuracy  of  +  0.1  percent)  and  a  minimum  conversion  rate  of  1000 
words  per  second.  The  outputs  shall  be  compatible  with  the  strip  chart  recorder 
inputs. 

B.3. 2  Strip  Recorder 

An  eight  channel  strip  recorder  must  be  provided  for  recording  analog 
signals.  The  inputs  to  the  recorder  are  to  be  compatible  with  the  analog  output 
channels  provided  so  that  digital  data  may  be  used  to  monitor  analog  signals 
during  on-line  processing.  The  strip  chart  recorder  will  be  a  Sanborn  Model 
7708A,  with  8801A  preamplifiers  or  equivalent. 

B.3. 3  Waveform  Display 

The  waveform  displays  will  be  used  to  collect  and  display  signals  taken 
from  any  seismometer  and/or  test  points  determined  by  subarray  control 
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commands.  All  seismometer  and  test  point  data  will  be  collected  and  stored  in 
a  processor  which  will,  upon  receipt  of  manual  requests,  format  selected  portions 
of  the  collected  data  into  messages  that  will  fit  the  requirements  of  a  waveform 
display. 

B.3.3.1  Display 

Eight  traces,  with  up  to  500  sample  points  per  trace,  will  be  displayed  on 
a  CRT.  The  abscissa  will  represent  time  and  the  ordinate  will  represent  linear 
signal  strength  (amplitude)  of  each  sample  point.  A  vertical  line  cursor  may  be 
superimposed  on  the  waveform  display.  Each  waveform  may  be  processed 
independently.  Provision  shall  be  made  for  vertical  magnification  of  the  wave¬ 
forms.  Display  flicker  will  be  kept  to  a  minimum  by  matching  its  repetition 
rate  to  a  specific  phosphor.  The  rate  should  be  not  less  than  30  Hz. 

The  phosphor  should  be  selected  so  that  smear,  under  normal  computer 
room  lighting  conditions,  will  not  be  apparent. 

B.3.4  Beam  Display 

The  Beam  Display  will  be  used  to  display  signals  which  have  been  taken 
from  the  LASA  array  seismometers,  formed  into  a  beam  (measure  of  seismic 
energy  level  at  a  specific  inverse  phase  velocity  location),  and  converted  into  a 
word  that  will  display  the  beam's  intensity  level.  A  large  number  of  beams  will 
be  formed  and  displayed  on  the  CRT.  Thus,  the  seismic  activity  at  a  large 
number  of  contiguous  inverse  phase  velocity  locations  will  be  displayed  at  the 
same  time. 

An  external  processor  will  collect  and  store  seismic  data,  form  the  beams, 
and  format  them  into  messages  that  will  fit  the  requirements  of  a  beam  display. 

B.3.4.1  Display 

A  32  by  32  matrix  of  beams  will  be  displayed  on  the  CRT.  Each  matrix  element 
will  be  the  position  of  a  beam  assigned  to  a  specific  inverse  velocity  space  location. 


Each  beam's  relative  intensity  (energy*)  level  will  be  indicated  by  the  spot 
intensity  (CRT  beam  current)  at  the  XY  coordinate  point  assigned  to  it.  Thus, 
light  output  at  a  beam's  coordinate  will  be  a  function  of  the  beam's  seismic 
energy  level.  The  functional  relationship  may  be  linear,  logarithmic,  or  floated 
logarithmic. 

Beam  size  must  be  adjustable  so  that  beams  do  not  overlap  and  at  the  same 
time  are  spaced  no  more  than  15  percent  of  the  beam  width  apart. 

Display  flicker  must  be  kept  to  a  minimum  by  matching  its  repetition  rate 
to  the  selected  phosphor.  The  rate  should  be  no  less  than  30  Hz. 

The  phosphor  should  be  selected  so  that  smear,  under  normal  computer 
room  lighting  conditions,  will  not  be  apparent. 

Update  rate  (fresh  data)  can  occur  at  a  maximum  rate  of  10  Hz. 

B.3.4.2  Manual  Inputs 

Manual  Inputs  (MI)  are  required  to  instruct  the  external  processor  on  what 
kinds  of  display  data  it  should  format  and  store  in  the  display  storage  element. 

An  MI  keyboard  overlay  providing  the  following  functions  is  required: 

a.  Select  the  beams  that  will  be  displayed. 

b.  Expand  the  display  around  some  specific  beam. 

c.  Freeze  a  selected  pattern  at  some  point  in  time. 

d.  Play  back  some  previous  pattern. 

e.  Select  the  functional  relationship  between  beam  seismic  energy  and 
display  brightness. 

f.  Select  and  display  a  specific  seismometer  output  at  some  specific 
trace  location. 

g.  Expand  a  selected  trace  a  specified  number  of  times. 

h.  Display  a  cursor  mark  on  a  selected  trace  or  beam;  shift  the  cursor 
in  conjunction  with  an  MI  entry  control,  relative  to  other  beams  or 
traces. 


B-5 


i.  Format  maintenance  instructions  that  can  be  used  to  control 
maintenance  activity  at  the  array. 

j.  Enter  a  seven  bit  pedestal  for  selecting  threshold  level. 

B.3.4.3  X-Y  Plotter  and  TWX  Facility 

The  equipment  must  include  an  X-Y  plotter  and  the  facility  to  operate  a 
remote  teletype  terminal. 

B.3.5  Array  Calibration  and  Maintenance  Facilities 
B.3  5.1  Calibration 

Calibration  of  the  seismometers  shall  be  accomplished,  a  subarray  at  a 
time,  by  transmitting  a  command  under  computer  control  to  the  subarray  which 
places  a  sine  wave  of  known  frequency  or  a  pseudo-random  waveform  of  known 
amplitude  on  each  sensor  calibration  coil.  The  data  received  from  each  seis¬ 
mometer  of  the  array  is  the  sum  of  the  known  test  waveform  and  any  ambient 
vibrations  occurring  at  that  time,  and  is  analyzed  by  the  computer  to  determine 
the  calibration  of  the  individual  seismometer  as  to  amplitude  and  frequency  of 
response. 

B.3.  5. 2  Maintenance 

a.  Manual — The  maintenance  console  shall  be  provided  with  eight  channels 
of  digital  to  analog  conversion,  eight  channels  of  strip  chart  recording,  and  a 
six  to  eight  trace  waveform  display  which  can  be  utilized  to  observe  and  record 
various  array  output  waveforms .  Under  manual  control,  by  the  setting  of 
switches,  any  subarray  shall  be  selectable  and  any  word  of  the  message  received 
from  that  subarray  shall  be  selectable  for  presentation  as  a  waveform  display 
or  made  available  at  any  one  of  the  eight  digital  to  analog  output  channels.  Each 
D/A  output  channel  may  be  connected  to  a  strip  chart  recorder  channel.  Keys 
shall  be  provided  to  designate  a  six  bit  command  to  initiate  its  transmission  to 

g 

a  selected  subarray  to  effect  the  various  tests  available  in  the  subarray  equipment. 
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b.  Computer  Controlled  Maintenance — By  diagnostic  programming,  any 
subarray  shall  be  selectable  to  receive  any  of  the  test  commands.  The  data 
received  from  the  subarray  may  be  recorded  on  the  strip  recorder  as  an  analog 
waveform,  may  be  displayed  on  a  waveform  scope,  or  an  analysis  may  be  made 
as  to  amplitude  and  frequency  by  the  computer  to  ascertain  the  validity  of  the 
seismometer  outputs.  The  significant  results  are  then  recorded  by  the  log-out 
printer. 

Under  program  control  the  American  Telephone  &  Telegraph  Co.  type 
X303A-10  5  modems  at  the  computer  will  be  connected  to  provide  a  loop  test 
capability,  utilizing  computer  generated  data  transmission  as  the  source  of  the 
messages. 

Additional  loop  test  capabilities  shall  be  provided  at  the  subarray  level, 
with  a  computer  generated  test  command  resulting  in  generation  of  subarray 
test  signals. 
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Appendix  C 


SEISMIC  DATA  COMPRESSION 


C.l  INTRODUCTION 

Present  LASA  instrumentation  consists  of  analog  sensors  with  corresponding 
amplifiers  capable  of  operating  throughout  a  72  dB  dynamic  range.  Precision 
analog  communications  are  feasible  for  short  distances,  but  digital  techniques 
are  more  favorable  when  traversing  long  distances  with  a  significant  composite 
data  rate.  This  has  led  to  a  subarray  concept,  where  an  ensemble  of  analog 
signals  from  local  sensors  are  collected  at  a  central  instrumentation  vault, 
band  limited,  sampled,  multiplexed,  digitized,  and  subsequently  transmitted  to 
the  LASA  communications  terminal.  A  fourteen  bit  fixed  point  conversion  is 
consistent  with  the  analog  instrumentation  performance  characteristics  and 
establishes  the  present  fundamental  data  rate  of  280  baud  per  instrument  when 
samples  are  acquired  at  a  20  Hz  rate.  Since  the  LASA  contains  525  seismom¬ 
eters,  a  147  kilobaud  (Kb)  communications  capability  is  required  to  acquire  di¬ 
rectly  the  short  period  instrument  data.  Secondary  instrumentation  and  overhead 
for  synchronization  and  error  control  has  extended  this  burden  to  201.6  Kb. 

This  is  presently  implemented  as  twenty -one  19.2  Kb  channels  which  impose  a 
403.2  Kb  facility  requirement. 

The  potential  advantages  of  processing  at  a  remote  location  creates  interest 
in  extending  the  LASA  communications  network.  Continued  utilization  of  the 
403.2  Kb  composite  channel  is  obviously  impractical,  and  even  the  base  data 
rate  of  147  Kb  constitutes  a  major  communication  cost.  A  survey  of  data  com¬ 
pression  techniques  is  appropriate  to  establish  the  feasibility  of  remote  process¬ 
ing  and  to  derive  more  efficient  data  acquisition  techniques  for  large  arrays. 
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C .  2  DELTA  MODULATION 


In  many  applications,  a  significant  data  reduction  has  been  achieved  by 
employing  delta  modulation;  i.e. ,  transmitting  the  signal  change  encountered 
during  the  sample  interval.  Characteristics  of  this  technique  are  presented  by 
Figure  C-l,  which  illustrates  the  word  length  reduction  potential  when  the  re¬ 
quired  sampling  rate  is  excessive  for  the  predominant  signal  components.  Note 
that  a  ratio  greater  than  6  samples  per  cycle  is  required  to  realize  a  range 
compression.  This  situation  is  often  encountered  in  broad  band  channels.  Prin¬ 
cipal  seismic  signal  spectrum  components  have  been  observed  at  1.2  Hz,  and 
are  conjectured  to  exist  as  high  as  1.4  Hz.  With  the  present  20  Hz  data  acquisi¬ 
tion,  14.3  samples  per  cycle  must  be  accommodated  which  permits  the  deletion 
of  only  one  bit  if  no  loss  in  performance  is  to  be  incurred.  This  7  percent  re¬ 
duction  is  significant,  and  would  presumably  justify  the  additional  instrumenta¬ 
tion  complexity  if  a  20  Hz  sampling  rate  was  essential  to  the  subsequent  process. 

C .  3  DATA  SAMPLING 

A  direct  method  of  data  compression  for  transmission  is  the  reduction  of 
the  sampling  rate  so  as  to  approach  the  theoretical  Nyquist  rate  requirement, 
followed  by  the  use  of  reconstruction  interpolation  at  the  receiving  end  to  re¬ 
produce  the  data  continuum.  This  technique  is  well  suited  for  data  acquisition 
from  LASA.  It  requires  that  only  minimal  changes  be  made  in  the  existing 
system  to  reduce  the  sampling  to  a  submultiple  rate  prior  to  transmission  and 
allows  selective  reconstruction  of  the  higher  sample  rate  data  within  a  general 
purpose  computer  when  additional  phase  resolution  is  required. 

If  the  sampling  rate  is  chosen  to  be  greater  than  twice  the  highest  data 
frequency,  sampling  theory  states  that  the  original  waveform  can  be  completely 
reconstructed.  Consider  the  signal  f(t)  with  the  Fourier  transformation: 
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Figure  C-i.  Delta  Modulation  Characteristics 
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such  that  the  signal  is  bandlimited  ( |  F(co)  |  =  0  for  |  oj  |  ^ 

Defining  the  sampling  frequency  as  2o^  and  the  time  interval  between  sam¬ 
ples  as  T,  the  sampling  frequency  would  then  be  chosen  such  that: 

Cbs  —  ^  Wi*  (C-2) 

Then  f(t)  can  theoretically  be  recovered  from  its  sampled  values  f(nT)  by  apply¬ 
ing  the  expression: 

®  sinco2(t-nT) 

f<‘>  -  n£,  f<nT>  Ml(t-nT)  <C-3> 

This  transform  corresponds  to  an  ideal  filter  that  passes  without  attenuation  all 
frequencies  below  tc2  and  completely  rejects  frequency  components  introduced 
by  the  sampling  process  which  exceeds  oj2. 

A  bandlimited  signal  cannot  be  realized,  and  some  energy  will  always  be 
encountered  beyond  the  capacity  of  the  sampled  data  channel  ( |  F (cc )  |  /■  0 
for  |  co|  s  co2) .  Even  if  the  sampling  rate  is  chosen  to  be  greater  than  twice  the 
highest  frequency  of  interest,  the  signal  cannot  be  perfectly  reconstructed,  for 
new  components  corresponding  to  the  data  spectrum  translated  by  the  harmonic 
structure  of  the  periodic  sampling  impulse  train  are  present  in  the  band  of 
interest.  This  aliasing  error  constitutes  the  theoretical  limit  on  data  recon¬ 
struction  fidelity.  However,  additional  error  is  imposed  by  the  practical  limi¬ 
tation  on  the  number  of  data  points  employed  by  the  reconstruction  algorithm, 
and  this  effect  is  examined  in  Appendix  C.4. 

In  the  present  LASA  system,  a  low -pass  filter  is  included  in  each  seis¬ 
mometer  channel  to  severely  attenuate  data  components  exceeding  5  Hz.  The 
amplitude  response  of  the  filter  and  the  alias  noise  introduced  into  the  data  for 
10  Hz  and  20  Hz  sampling  rates  are  shown  in  Figure  C-2  where  the  input  has 
been  modeled  as  wideband,  uniform  spectral  density  data.  Lowering  the  cutoff 
frequency  of  the  filter  reduces  the  amount  of  noise  aliased,  or  folded  into  the 
data  spectrum.  If  the  aliasing  error  at  a  reduced  sampling  rate  is  deemed  un¬ 
acceptable,  alteration  of  the  present  filters  should  be  considered. 
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Fi rnn-Q  £-2.  Frequency  Plane  Aliasing 


A  useful  technique  for  computing  the  upper  bound  on  the  aliasing  error  has 
been  developed  by  Papoulis8  and  modified  by  Stickler9.  The  maximum  error 
|e(t)  |  in  a  reconstructed  function  f(t)  is  bounded  by  the  inequality: 

/.OD 

|e(t)  |  <  |  /  | F(co)  |  dco .  (C-4) 

J(jJz 

If  F(co)  is  bandlimited  (F(co)  =  0  for  |co  |  s  co0>  002),  the  expression  becomes: 

?  fW  o 

|e(t)  |  <  -  I  |F(W)|da).  (C-5) 

J  coz 

This  bound  can  be  used  as  a  guide  for  estimating  the  effect  of  sample  rate  re¬ 
duction. 

In  the  following  analysis,  data  with  a  20  Hz  sample  rate  is  generated  from 
the  10  Hz  data.  The  original  20  Hz  data  is  utilized  as  a  reference  to  evaluate 
the  performance  of  several  reconstruction  algorithms.  Thus,  the  input  can  be 
considered  bandlimited  at  10  Hz  which  corresponds  to  the  capacity  of  the  refer¬ 
ence  data.  An  upper  bound  on  the  aliasing  error  can  be  expressed  as: 


(C-6) 


Observations  from  seismometer  AO  of  subarray  AO  were  employed  for  this 
analysis.  The  first  record  consisted  of  51  seconds  of  seismic  noise  and  the 
second  was  a  25  second  period  of  data  recorded  for  the  Kamchatka  Earthquake 
of  April  8,  1966.  The  following  upper  bound  of  aliasing  errors  have  been  cal¬ 
culated  from  spectral  density  estimation  data: 


DATA 


ESTIMATED  ALIASING  ERROR  BOUND 


Noise  Record 


-23.7  dB  relative  to  the  peak  data  magnitude 


Kamchatka  Earthquake  -31.9  dB  relative  to  the  peak  data  magnitude 
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C . 4  DATA  RECONSTRUCTION 


Resolution  requirements  are  often  misconstrued  as  an  essential  sampling 
constraint,  but  actually  they  constitute  a  blend  of  sampling  and  reconstruction 
requirements.  When  communications  represent  a  major  system  consideration, 
reconstruction  techniques  should  be  considered  for  acquiring  the  desired  record 
resolution.  The  rudimentary  theoretical  assumptions  are  only  approximated  in 
practice  and,  since  the  most  uncertain  data  points  lie  midway  between  the  avail¬ 
able  datums,  candidate  interpolation  algorithms  for  estimating  these  points  have 
been  experimentally  evaluated. 


Three  methods  of  doubling  data  density  by  using  linear  symmetrical  com¬ 
binations  of  up  to  six  consecutive  points  on  each  side  of  the  estimated  point  were 
simulated.  The  base  10  Hz  data  was  generated  by  selecting  every  other  value 
of  20  Hz  data,  and  hence  the  interpolated  values  can  be  compared  to  the  orig¬ 
inally  observed  values.  The  three  methods  are  summarized  as  follows: 


a.  Lagrangian  Polynomial  Interpolation 

Given  a  set  of  m+1  equidistant  data  values  f(KT),  K=j, . . . ,  j-m,  there 
is  only  one  polynomial  0(t)  of  degree  m  such  that  0(KT)  =  f(KT)  for 
j-m  £  K  ^  j.  0(t)  can  be  generated  from  the  Lagrangian  polynomials10. 

Polynomials  of  degree  1  and  3  were  evaluated  for  reconstruction  inter¬ 
polation.  The  degree  1  or  linear  interpolation  algorithm  for  the  point 
equidistant  between  two  samples  is  given  by: 

X(nT)  -  X|  (n+1)T1  +  X[  (n-l)T]  (C.7) 

2 


The  degree  3  interpolation  algorithm  for  the  point  in  the  center  of  the 
four  samples  is  given  by 

£(nT)  =  9{X[(n-l)T]  +X[(n+1)T^  -  jx[(n-3)T]  +X[(n+3)T][  _  (C_g. 

b.  Sinx/x  Reconstruction  Interpolation 

The  ideal  sinx/x  reconstruction  is  approximated  by  using  N  points  on 
each  side  of  a  sample  interval  to  reconstruct  the  midpoint.  The 
algorithm  for  this  case  is  given  by: 


_  S  2  {x[(n-2k+l)T]  +  X[(n+2k-l)T]  |  <-l)k+1 

X(nT)  '  I,  - (2^1)i 


where  N  =  1,  2,  3,  4,  5,  6. 


(C-9) 
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c.  Optimum  Discrete  Filter  Interpolation 

The  optimum  or  Wiener  interpolation  filter  used  here  minimizes  the 
interpolation  error  in  the  mean  square  sense  under  the  constraint  that 
the  interpolated  values  are  a  linear  combination  of  the  sample  values: 

a  N-l 

X(t)  =  £  h,  X'(n+k)T]  for  the  interval  (n-l)T  <  t  <  nT,  (C-10) 

k=-N 

A 

where  X(t)  is  the  interpolated  value  of  X(t)  at  time  t.  The  mean  square 
error  function  is 

MSE  =  E  j  [X(t)-X(t)p  J  ,  (C-ll) 

A 

where  E  denotes  the  expected  value.  Substituting  for  X(t)  and  mini¬ 
mizing  the  resulting  expression  with  respect  to  the  coefficients  h^ 
yields,  under  the  assumption  that  X(t)  is  a  wide  sense  stationary  ran¬ 
dom  process  with  autocorrelation  Rxx(t),  the  following  set  of  equations 
which  specify  the  proper  coefficient  values: 


N-l 

£ 

k=-N 


VW<'-k>Tl 


=  R 


XX 


UT-t) 


(C-12) 


and  l  =  -N,  -N+l, . . . , N-2.N-1. 

These  equations  were  evaluated  for  values  of  N  from  one  through  six 
using  flat  power  spectrum  data  models  from  zero  to  3  Hz  and  from  0.  5 
to  2.5  Hz. 

This  technique  could  be  extended  to  yield  the  best  interpolation 
estimate  (in  the  mean  square  sense)  of  the  signal  for  the  case  when  sig¬ 
nal  and  noise  autocorrelation  functions  are  individually  described. 
Assuming  that  the  signal  and  noise  are  uncorrelated,  the  coefficients 
hfc  are  specified  by  the  following  set  of  equations: 


2hkiRss[^-K>Tl  +  Rnn[^_K>T]  1  =  Rss^T't>  (C-13) 

and  t  =  -N,  -N+l, . . .  ,N-2,N-1  where  Rss(t)  and  Rnn(T)  are  the  auto¬ 
correlation  functions  of  the  signal  and  of  the  noise  respectively. 


The  optimum  interpolation  error  for  the  0.5-2. 5  Hz  spectrum  case  was 
within  approximately  1  dB  of  the  error  for  the  0-3  Hz  spectrum  case  for  the 
same  algorithm  data  base.  Only  the  latter  is  graphically  represented.  Figures 
C-3  and  C-4  show  the  ratio  of  maximum  error  magnitude  to  maximum  data  mag¬ 
nitude  when  each  of  the  algorithms  is  applied  to  the  Kamchatka  earthquake  and 
seismic  noise  records.  Figures  C-5  and  C-6  show  the  ratio  of  the  power  of  the 
error  to  the  total  data  power  in  dB  for  the  two  records.  The  power  density  for 
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Kamchatka  Earthquake  (4-8-66)  As  Indicated 
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Figure  C-3.  Reconstructed  Signal  Error  Bound 


20  Hz  Data  Estimated  From  10  Hz  Data 
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Figure  C-4.  Reconstructed  Noise  Error  Bound 


Kamchatka  Earthquake  (4-8-66)  As  Indicated 
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Figure  r-5.  Reconstructed  Signal  Error  Power 


20  Hz  Data  Estimated  From  10  Hz  Data 
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Figure  C-6.  Reconstructed  Noise  Error  Power 


the  error  between  the  actual  values  X(nT)  and  the  estimated  values  X(nT)  gen¬ 
erated  by  the  3rd  order  polynomial  algorithm  are  illustrated  by  Figures  C-7 
and  C-8. 

Observation  of  these  experimental  results  reveals  a  reconstruction  per¬ 
formance  which  introduces  an  insignificant  additional  noise  component  consistent 
with  the  anticipated  aliasing  bound.  Thus,  an  excessive  sampling  rate  is  pres¬ 
ently  implemented,  and  the  compatible  submultiple  of  10  Hz  appears  sufficient  to 
retain  the  message  fidelity.  This  50  percent  communication  requirement  re¬ 
duction  precludes  the  desirability  of  delta  modulation,  for  the  principal  compo¬ 
nent  sampling  ratio  is  no  longer  adequate  to  provide  a  significant  range 
compression. 

C.  5  DATA  ENCODING 

Additional  data  compression  can  be  achieved  by  examining  the  trade-offs 
between  dynamic  range  and  signal  resolution  for  a  floating  point  data  word 
instead  of  the  present  fixed  point  format.  In  the  fixed  point  format,  the  resolu¬ 
tion  of  the  signal  is  a  function  of  its  magnitude;  large  magnitude  signals  have 
excessive  resolution  to  achieve  a  large  dynamic  range.  By  scaling  the  signal 
to  operate  near  the  maximum  value  acceptable,  a  minimum  resolution  over  the 
entire  dynamic  range  can  be  achieved  with  fewer  data  bits .  Implementation 
complexity  can  be  minimized  by  retaining  the  two's  complement  notation  for 
negative  numerical  values  and  utilizing  integer  powers  of  two  as  the  arithmetic 
base. 

Figures  C-9,  C-10,  and  C-ll  show  the  relationship  between  the  resolution 
band  and  dynamic  range  as  a  function  of  word  size  for  base  2,  base  4,  and 
base  8  floating  points.  The  word  size  includes  the  sign  and  the  exponent  fields. 
Figure  C-12  shows  the  relationship  of  the  minimum  resolution  and  dynamic 
range  for  the  three  formats.  Base  2  and  base  4  have  the  same  minimum  res¬ 
olution  for  a  given  word  length  and  dynamic  range,  but  base  4  has  a  higher 
maximum  resolution.  In  the  dynamic  ranges  and  word  lengths  of  interest, 
base  2  and  base  4  are  both  better  than  base  8. 
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Figure  C-7.  Signal  Error  Spectrum-3rd  Order  Polynomial 
Reconstruction  Algorithm 
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Frequency  (Hz) 

Figure  C-8.  Noise  Error  Spectrum— 3rd  Order  Polynomial 
Reconstruction  Algorithm 
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Figure  C-9.  Floating  Point  Base  2  -  Encoding  Characteristics 


NUMBER  OF  BITS 
(INCLUDING  SIGN 
AND  EXPONENT 
FIELDS) 


1 


I - i  i - 1 - 1  . .  i - 1 

cn  o  co  >o  ^  <nj  ° 

(gp)  ONvg  Nonmosaa 


C-17 


Figure  C-10.  Floating  Point  Base  4  -  Encoding  Characteristics 
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Figure  C-ll.  Floating  Point  Base  8  -  Encoding  Characteristics 
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Dyncamic  Range  (dB) 

Figure  <^-12.  Relative  Encoding  Characteristics 


The  most  suitable  selection  for  this  specific  application  is  base  4  modified 
floating  point  encoding.  A  10-bit  word  is  capable  of  representing  a  bi-polar 
value  throughout  a  78  dB  dynamic  range  with  a  large  signal  resolution  of  0.27  dB. 
The  introduced  quantization  noise  is  at  least  30  dB  below  the  largest  power  of  in¬ 
stantaneous  signal  or  background  noise,  and  192  numerically  illegal  codes  are 
available  for  fault  signaling  purposes.  This  technique  is  directly  compatible 
with  the  present  LASA  instrumentation  dynamic  range,  but  no  requirement  for 
a  quantum  scaled  40  dB  below  the  inherent  noise  is  recognized.  Since  the 
processing  techniques  which  attenuate  the  background  noise  are  also  effective  on 
quantum  noise,  a  rescaling  by  approximately  20  dB  should  be  considered  to  ex¬ 
tend  the  system's  linear  range  to  magnitude  seven  events.  This  recommenda¬ 
tion  is  illustrated  by  Figure  C-13,  which  also  depicts  the  effective  array  reso¬ 
lution  extension  which  prevails  as  a  product  of  the  beam  forming  process.  Less 
than  0.01  dB  of  additional  uncorrelated  background  noise  is  incurred  at  low 
levels,  and  even  more  extreme  rescaling  may  be  advisable  when  the  subsequent 
processes  are  firmly  established. 

This  communication  concept  extends  the  system  performance  while  reducing 
the  data  rate  to  100  baud  per  short  period  instrument.  Signal  analysis  has 
demonstrated  that  the  inner  instrument  ring  of  each  LASA  subarray  can  be  de¬ 
leted  without  adverse  effects.  Primary  seismic  data  is  then  reduced  to  39.  9 
Kb,  and  can  be  easily  transmitted  through  the  relatively  economical  Bell 
TELPIPE  with  an  adequate  22  percent  margin  for  secondary  instrumentation 
and  overhead.  These  techniques  also  illustrate  the  feasibility  of  employing 
voice  grade  lines  for  the  data  acquisition  network  coupling  the  subarrays. 
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Figure  C-13.  Seismometer  Digital  Communications  Scaling 


Appendix  D 

DETECTION  PROCESSOR  SCALING 

D.  1  INTRODUCTION 

Because  of  the  finite  word  length  limit  in  a  digital  computer,  there  are 
unavoidable  errors  due  to  quantization  of  the  input,  of  the  coefficients  of  the 
digital  filters,  and  of  the  results  of  computation.  A  study  of  these  quantization 
(or  truncation)  errors  is  important  to  ascertain  the  adequacy  of  implementation, 
to  permit  a  prudent  assignment  of  variable  scaling,  and  to  evaluate  the  numer¬ 
ical  error  present  at  each  process  node. 

D.  2  SCALING 

The  effects  of  quantization  are  modeled  as  additive  random  noise.  These 
contributions  are  combined  with  the  estimated  data  noise  encountered  at  each 
stage  of  the  Detection  Processor,  and  the  composite  noise  propagation  is  traced 
through  the  system.  The  results  are  in  reasonable  agreement  with  preliminary 
experimental  observations. 

D.  2. 1  Seismometer  Data 

The  output  of  a  single  seismometer  is  digitized  as  a  13-bit  word  plus  sign 
at  a  quantum  level  (Q)  of  0.  028  nm.  The  quantization  error  is  assumed  to  be  a 
random  variable  which  is  uniformly  distributed  and  uncorrelated  from  sample 
to  sample.  Since  two's  complement  representation  is  used  within  the  computer, 
the  quantization  error  bias  (B)  is  and  the  variance  is  Q^/12. 

The  natural  background  noise  is  observed  to  be  a  variable  with  zero  mean 
and  a  standard  deviation  a  of  approximately  3  nm,  or  roughly  100  times  the 
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quantization  level.  The  noise  power  due  to  quantization  is  therefore  five  orders 
of  magnitude  lower  than  the  power  due  to  the  background  noise,  and  the  former 
may  be  neglected  for  further  practical  consideration. 


Summarizing,  at  the  seismometer  level, 

Q  =  2.  8  x  10-2  nm 
B  =  -Q/2  =  -1.4  x  10-2  ^ 


(D-l) 


a  =  3  nm  =  1.07  x  102  Q. 

D.  2.  2  Subarray  Beamforming 

Although  the  microcode  for  this  process  has  the  capability  to  scale  down 
large  magnitude  inputs,  it  is  not  considered  in  this  model.  A  prescaling  right 
shift  of  five  or  fewer  binary  digits  will  not  effect  the  process  parameter  selec¬ 
tion.  The  subarray  beamforming  is  carried  out  by  properly  delaying  and  sum¬ 
ming  the  outputs  of  the  25  seismometers  in  the  subarray.  Since  addition  only 
is  involved  in  this  implementation,  no  other  computational  error  is  introduced. 

At  the  subarray  beam  level,  there  is  a  gain  in  signal-to-noise  ratio.  For 
Longshot,  the  observed  unfiltered  subarray  gain  was  4.5  dB,  or  a  factor  of  1.68. 
This  broad  band  apparent  gain  is  utilized  to  estimate  the  residual  noise  at  the 
conclusion  of  the  process.  Thus,  at  the  subarray  beam  level, 


2.8  x  1(T2 


1. 12  x  10"3  nm 


B  =  -1.4  x  10-2  nm  =  -1.  25  x  10  Q 


(D-2) 


D.  2.3  Recursive  Filtering 

The  filter  selected  for  initial  detection  processor  evaluation  is  a  third 
order  bandpass  Butterworth  formulation  characterized  by  the  parameters 
appearing  in  Table  D-l. 
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TABLE  D-l 

FILTER  CHARACTERISTICS 


Parameter 

Value 

Sample  Period 

T 

0.1  Seconds 

Lower  Corner 

Frequency 

h 

0.9  Hz 

Upper  Corner 

Frequency 

h 

1.4  Hz 

Output/Input 

Bound  Ratio 

5>i 

1.59 

n=0 

Implementation 

Bias  Magnification 

Jd„ 

12.38 

n=0 

Implementation 

Noise  Magnification 

3264.47 

n=0 

The  filter  function, 


H(Z)  -  - -  K(l-Z^)3 -  ,  (D-3) 

1+bjZ  +  b2  Z  ^  +  .  .  .  +  bg  Z  6 

is  obtained  by  implementing  the  difference  equation  in  accordance  to  Figure  D-l 
Theoretical  values  of  the  coefficients  are  listed  in  the  first  column  of  Table  D-2 
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Figure  D-l.  Recursive  Filter 


TABLE  D-2 

FILTER  COEFFICIENTS 


Parameter 

Value 

Theoretical 

Machine  Scaled 

Implemented 

K 

0.002898 

3 

bl 

-4.081417 

-8358.74 

-8359 

b2 

7.955923 

16293.73 

16294 

b3 

-9.092820 

-18622.09 

-18622 

b4 

6.444996 

13199.35 

13199 

b5 

-2.677543 

-5483.61 

-5484 

b6 

0.532075 

1089.69 

1090 

The  scaling  factor  of  2-^  relating  the  theoretical  and  machine  scaled  values 
is  implied  by  the  objective  to  minimize  the  encountered  quantization  while 
accommodating  the  largest  denominator  coefficient  (bg  )  within  the  15  bit 
fixed  point  computation  word  Integer  numbers  must  be  employed  to  approxi¬ 
mate  the  machine  scaled  values. 

The  effects  of  fixed  point  word  saturation  and  the  resulting  phenomenon  of 
limit  cycles  have  been  experienced  during  processor  development.  This  may 
be  prevented  by  sufficiently  restricting  the  input  so  as  to  preclude  an  excessive 
output  value.  The  output  to  input  bound  ratio  corresponds  to  the  sum  of  the 
magnitudes  of  the  impulse  response  of  the  filter.  In  this  instance, 

OO 

is  1.59,  thereby  establishing  an  input  magnitude  bound  of 

n=0 


1.25  x  2^4  quantum  units  if  linear  operation  is  to  be  assured.  The  filter 
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scaling  coefficient  (K)  restriction  is 


K  < 


0.002898  x  211 

1.59 


3.73 


which  utilizes  input  saturation  to  preclude  saturation  of  the  filter  process.  The 
largest  permissable  integer  is  3,  and  thus  the  transfer  function  numerator  is 
implemented  as: 

-11  9  3 

A(z)  =  3  x  2  11  (1  -  Z"2)  . 

This  establishes  a  new  quantization  level  (within  the  filter  pass  band): 


Q 


0.00112  x  0.002898  „  „„  ,„-3 

- — -  =  2.22  x  10  nm 

3  x  2”11 


(D-4a) 


Figure  D-2  is  a  plot  of  the  filter  response  calculated  with  the  theoretical 
coefficients.  The  response  calculated  with  the  implemented  coefficients 
(Figure  D-3)  gives  an  almost  identical  curve  that  has  been  shifted  downward  as 
a  result  of  scaling  the  numerator  to  prevent  filter  saturation.  The  normalized 
difference  between  the  response  curves  is  illustrated  in  Figure  D-4. 

Truncation  errors  encountered  in  the  filter  process  contribute  to  the  output 
noise  level.  Consider  shifting  a  word  of  m  digits  right  r  places  and  then  trun¬ 
cating  to  create  an  effective  word  length  of  m-r  digits.  The  truncation  error 
due  to  the  right  shift  is  equivalent  to  the  difference  between  the  prevalent  quan¬ 
tization  errors  associated  with  the  m  bit  and  m-r  bit  word  lengths.  The 
bias  and  variance  of  the  error  component  are  related  to  the  quantum  level: 

Br  =  -(l-2“r)  Q/2 

or2  =  (l-2“2r)  Q2/l2, 

where  Q  represents  the  quantization  level  after  the  right  shift.  The  second  term 
in  each  of  the  above  expressions  is  negligible  when  evaluating  the  effects  of  the 
filter  scaling  shift  (r=  11 ) . 
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Figure  D-2.  Frequency  Response  for 
Theoretical  Filter 
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Figure  D-3.  Frequency  Response  for 
Actual  Filter 
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Figure  D-4.  Deviation  of  Frequency  Response 
for  Actual  Filter  from  Theoretical 
Filter 
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Statistical  characteristics  of  the  truncation  error  generated  in  the  filter 
recursive  loop  are: 


Bf  =  (-Q/2)  ^  dn  =  -6.19  Q 

n=0 


(Q2/12 )  7  dn2  =  1.65  x  10  Q 


n=0 


where  dn  is  the  unit  impulse  response  of  the  filter  denominator  implementation. 
The  errors  are  assumed  to  be  uniformly  distributed  and  uncorrelated  from 
sample  to  sample. 

No  bias  component  is  propagated  through  the  filter  since  a  null  transfer 
relation  exists  at  zero  frequency,  as  is  verified  by  the  final  value  theorem: 

OO 

y  cn  -  lim  H(Z)  =  0. 

Lt  z->  l 

n=0 


For  Longshot,  the  observed  gain  in  the  signal-to-noise  ratio  in  the  recursive 
filter  is  approximately  18.5  dB,  or  an  amplitude  ratio  of  8.35.  To  establish  an 
output  noise  estimate,  the  signal  is  assumed  to  be  coherent  and  the  gain  due 
entirely  to  noise  reduction.  The  standard  deviation  of  this  propagated  output 
noise  component  is: 

(1.59*  103  Q)  (1.12x10-3_)  „  9.62xl0Q. 

8.35  (2.  22  x  10“3) 


The  composite  noise  output  estimate  can  thus  be  expressed  as: 

B  =  -6.19  Q 

<7  -  {(1.65  x  10  )2  +  (9.62  x  10)2  }1/  2  Q  =  9.8  x  10  Q  (D-4b) 
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D.2.4  Filter  Output  Scaling 

The  output  of  the  recursive  filter  is  right  shifted  four  bit  positions  so  as  to 
best  match  the  dynamic  range  of  the  LASA  beamformer  and  the  preceding 
process.  This  alters  the  quantum  level: 

Q  =  (2.22  xlO"3)  x  24  =  3.55  x  10-2  nm.  (D-5a) 

Additional  error  is  introduced  by  the  shift,  and  the  bias  and  standard  deviation 
are  now: 


B 


-  6,^Q  -f  (1-2-4)  =  -8.56  x  IQ"1  Q 


I  (9.8  x  10  ) Q 


16 


^2  (i"2"8)  =  6.11  Q. 


(D-5b) 


D.2.5  LASA  Beamforming 

The  LASA  beams  are  formed  from  the  delayed  sum  of  21  subarray  beams. 
Again,  there  are  only  additions  involved.  The  gain  in  LASA  beamforming  is 
approximately  11  dB,  or  an  amplitude  factor  of  3.55. 

Q  =  3.55  x  10"2/21  =  1.69  x  l(f3  nm 

B  =  -(8.56  x  10-1)  x  21  Q  =  -1.798  x  10  Q  (D-6) 


Q  =  6,11  x  21  Q 
3.55 


3.61  x 10  Q 


The  predictable  bias  should  be  removed  prior  to  rectification  so  as  to  mini¬ 
mize  the  post  detection  noise.  This  can  be  accomplished  by  implementing  a 
22nd  input  channel  which  inserts  a  compensating  positive  signal  of  18  quantum 
units.  The  estimated  residue  is  only  2  x  10-2  Q,  or  about  3  x  10”^  nm. 
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D.2.6  Rectification 


A  sum  of  absolute  magnitudes  is  obtained  for  every  mutually  exclusive 
group  of  three  successive  samples,  and  this  value  is  right  shifted  two  bit  posi¬ 
tions.  The  subsequent  quantum  value  is: 

Q  =  |  (1.69  x  10-3)  =  2.25  xl0'3nm.  (D-7a) 


Assuming  a  Gaussian  distribution,  the  bias  of  the  output  of  a  full-wave 

nr 

linear  rectifier  with  unity  gain  is  evaluated  as  \j  —  times  the  standard  devia- 

7T 

tion  of  the  input  when  the  mean  is  zero11. 


B 


3 

4 


3.62  x  10  Q 


=  2.17  x  10  Q. 


(D-7b) 


The  standard  deviation  relationship  is  also  evaluated  as: 
(0.603)  (3.62  x  10)  Q  =  1.64  x  10  Q. 


(D-7c) 


D.2.7  Short  Term  Average 

An  accumulation  of  three  of  the  above  groups  corresponds  to  a  1.8  second 
short  term  window,  which  is  consistent  with  the  recommendations  of  the  post¬ 
detection  integration  analysis1.  This  establishes  the  final  output  character¬ 
istics  as: 

Q  =  |  (2.25x10-3)  =  7.5  x  10'4  nm 

B  =  3  (2.  17x10  )Q  =  6.51x10  Q  (D-8) 

o’  =  3  (1. 64  x  10 )  Q  =  4.  92  x  10  Q. 


D.2.8  Long  Term  Average 

The  background  noise  level  is  estimated  so  as  to  permit  the  application  of 
a  floating  threshold  criterion.  An  80-second  history  has  been  found  necessary 
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to  obtain  consistent  long  term  averages.  An  exponential  filter  of  the  form 


(D-9) 


has  been  implemented  to  preclude  the  requirement  for  maintaining  such  an 
extensive  history  within  the  processor.  This  is  accomplished  with  the  configu¬ 
ration  shown  in  Figure  D-5. 


The  factor  (1-2  )  represents  the  time  constant  of  the  exponential  filter. 
Figure  D-6  depicts  the  values  of  (3  versus  time  constant  normalized  to  the 
sample  period.  (3  was  chosen  to  be  6,  which  corresponds  to  a  time  constant  of 
113  seconds.  The  value  of  1  was  selected  for  a  so  that  the  bias  error  generated 
in  the  filter  would  alter  the  noise  estimate  by  only  0.1  dB. 

With  this  parameter  set,  the  dc  gain  of  the  filter  can  be  shown  to  be  32;  and 
the  quantization  is : 


-4 


7.5  x  10 
32 


2.34  x  10  0  nm. 


(D-lOa) 


Q 


Both  internally  generated  noise  and  propagated  noise  are  evident  at  the  output. 
Their  contributions  are  evaluated  as: 
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Figure  D-5.  Exponential  Filter 
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where: 
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Thus,  B  =  2. 05  x  103  Q 
and  <j  =  1.  39  x  103  Q 
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D.  3  SUMMARY 

These  quantizations  and  errors  for  the  Detection  Processor  implementation 
are  summarized  in  Figure  D-7.  The  standard  deviations  tabulated  are  based 
upon  the  noise  measurements  and  gain  performance  measured  for  Longshot. 
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Figure  D-7.  Detection  Processor  Scaling 
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EXPERIMENTAL  DISPLAY 

To  provide  rapid  recording  of  waveforms,  for  support  of  studies  during  the 
experimental  phases  of  the  LASA  Signal  Processing  system,  a  strip  chart 
recorder  interface  is  being  added  to  the  experimental  display  control  unit. 
Simultaneous  recording  of  multiple  time-varying  states  for  both  detection  and 
event  processing  during  on-line  operation  will  be  possible. 

The  strip  chart  recorder  will  utilize  the  experimental  display  interface  with 
the  System/360  through  the  2701  parallel  data  adapter.  (See  Figure  E-l.)  A 
maximum  of  eight  output  channels  will  be  available,  and  data  recording  will  be 
performed  at  a  faster  than  real-time  rate.  In  addition,  start  and  stop  of  the 
paper  transport  will  be  under  control  of  the  CPU.  Waveform  identification, 
starting  time,  scaling,  and  time  interval  can  be  provided  prior  to  each  record 
with  a  binary  coding.  The  width  of  the  time  interval  pulse  represents  a  known 
time  record  while  the  amplitude  represents  a  known  input  magnitude.  By 
shifting  the  magnitude  of  the  time  interval  pulse  the  same  as  data,  a  known 
amplitude  reference  is  provided.  An  example  of  a  typical  recording  is  shown 
in  Figure  E-2. 

A  data  buffer  consisting  of  the  following  information  has  been  defined: 
display  indicator  half-word,  channel  identification  sequence,  starting  time 
sequence,  time  interval  pulse,  followed  by  waveform  intensities  arranged 
sequentially  for  each  time  sample.  This  buffer  sequence  is  shown  in  Figure  E-3. 
The  high  order  byte  of  the  indicator  half-word  (bits  0  thru  7)  is  used  to  start 
the  strip  chart  recorder  paper  transport  and  indicate  the  number  of  channels 
being  recorded.  This  information  is  used  to  control  the  data  timing.  Once  the 
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Figure  E-l.  Experimental  Display 
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Figure  E-2.  Illustrative  Strip  Chart  Record 
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data  transfer  has  been  started,  the  operation  continues  as  a  background  program 
with  data  transfer  under  control  of  the  interface  timing.  The  number  of  half¬ 
words  indicated  in  the  indicator  byte  is  transferred  in  a  burst  at  a  rate  of 
approximately  200K  bytes  per  second.  After  transmission  of  the  last  half¬ 
word  in  the  burst,  the  2701  control  unit  will  not  respond  that  it  has  received 
the  data  until  it  is  time  to  transmit  a  new  burst  of  half-words.  The  time  between 
transmitting  bursts  must  be  sufficient  to  permit  the  strip  chart  recorder  writing 
pen  to  respond  to  the  change.  The  operator  will  be  able  to  control  the  time 
between  the  start  of  the  bursts  in  discrete  steps  from  8.33  milliseconds  to  0.2 
seconds.  For  records  with  low  sample  rates,  recording  time  will  be  much 
faster  than  real  data  time.  As  an  example,  a  two  minute  record  with  a  sample 
rate  of  five  samples  per  second  and  a  recorder  setting  for  an  interval  between 
the  start  of  data  bursts  of  25  milliseconds  would  require  15  seconds  to  record. 
The  record  time  is  independent  of  the  number  of  waveforms  recorded 
(1  through  8). 

Initiation  of  a  strip  chart  recording  operation  requires  the  presence  of  a 
data  buffer,  selection  of  the  strip  chart  recorder  mode,  setting  the  scaling 
control  and  activating  the  interrupt  pushbutton  on  the  display  control  panel. 

An  ATTENTION  is  presented  to  the  CPU,  which  shall  respond  by  scanning  the 
display  control  settings.  When  scaling  is  required,  the  data  in  the  data  table 
(all  except  the  indicator  half-word)  is  shifted  left  the  required  number  of  bit 
positions  while  retaining  the  proper  sign  for  the  data.  If  the  scaling  results 
in  overflow  of  data  bits,  the  magnitude  is  set  to  its  maximum  value  while 
retaining  its  original  sign.  This  shifted  data  is  used  to  assemble  a  new  data 
table  for  recording.  The  original  data  table  is  retained  in  case  the  operator 
requires  a  different  scaling  setting  for  data  that  has  overflowed  in  the  scaling 
operation.  After  completion  of  the  scaling  operation,  the  CPU  initiates  a 
WRITE  operation.  If  the  operator  changes  a  control  setting  during  the  write 
operation,  the  CPU  will  receive  a  DEVICE  END  and  ATTENTION  status  from 
the  2701  control  unit.  The  CPU  initiates  a  READ  operation  for  scanning  the 
control  settings,  performs  the  necessary  scaling,  and  restarts  the  data  transfer 
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from  the  beginning  of  the  new  data  buffer.  The  display  control  operation  has 
been  previously  described  in  conjunction  with  the  beam  display  operation5. 

The  selection  of  the  waveforms  being  displayed,  starting  time  for  the  data 
buffer,  and  generation  of  the  original  data  buffer  are  under  control  of  another 
I/O  device.  All  seismometers,  subarray  beams,  LASA  beams,  and  other 
desired  waveforms  will  be  assigned  an  identification  number  to  be  used  for 
selection  and  record  identification.  If  data  being  presented  has  been  recorded 
at  different  sample  rates  on  the  data  tape,  the  data  buffer  is  generated  to 
correspond  to  the  highest  sample  rate  with  either  the  value  at  the  last  sample 
time  or  an  interpolated  value  inserted  for  time  periods  where  data  does  not 
exist. 
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ARRAY -PRO CESSING  DESIGN  SYNTHESIS 

F.  1  INTRODUCTION 

The  purpose  of  this  appendix  is  to  identify  an  array-processing  design 
synthesis  procedure  upon  which  an  experimental  data  measurement  plan  can  be 
based.  Using  the  limited  data  base  processed  at  this  time,  certain  key  system 
design  parameters  can  be  measured  and  used  to  determine  expected  system 
performance,  b 2  ?  ,4 

It  is  intended  to  extend  the  synthesis  procedure  presented  here,  unify  the 
approach,  and  provide  substantially  more  experimental  verification.  In  addition, 
a  cost  analysis  along  the  lines  suggested  herein  is  under  study  and  may  be 
integrated  to  relate  cost  trade-offs  to  system  performance.  In  this  manner,  an 
array-processing  synthesis  technique  can  be  developed  pertinent  to  overall 
system  objectives,  optimally  configured  to  take  advantage  of  cost  trade-off's, 
and  expressed  in  terms  of  measurable  variables  to  enable  experimental 
verification. 

F.  2  THEORETICAL  CONSIDERATIONS 

The  communications  requirements,  system  gain,  and  location  accuracy  are 
of  interest  to  large  seismic  array -processing  design  synthesis. 

F.2.1  Data  Transmission 

The  dynamic  range  of  field  instruments  seldom  exceeds  50  dB,  and  the 
ratio  of  resolution  to  dynamic  range  is  usually  about  1  percent.  In  the  case  at 
hand,  we  are  interested  in  a  dynamic  range  of  well  over  70  dB  and  a  resolution 


F-l 


of  the  order  of  1/4  dB.  Since  the  instruments  will,  in  all  probability,  be  remotely 
located,  a  data  transmission  system  will  be  required.  Generally  speaking,  it 
is  desirable  to  reduce  the  data  rate  (bits  per  second)  to  a  minimum  value  while 
maximizing  channel  loading  within  available  capacity. 

The  dynamic  range  of  a  floating  point  modulo  4  communication  channel  and 
the  word  length  requirement  for  several  values  of  resolution  are  discussed  in 
Appendix  C.  This  data  formatting  technique  is  easily  coded  and  decoded  and  can 
be  supported  over  standard  voice  grade  telephone  lines.  It  is  noted  from  Figure 
F-l  that  a  10  bit  data  word  length  will  yield  a  dynamic  range  of  about  78  dB  with 
a  resolution  of  approximately  l/4  dB.  From  the  spectral  density  data  shown  in 
Section  F.3.3,  it  is  reasonable  to  suggest  a  sampling  rate  of  10  Hz,  which  is 
over  three  times  the  frequency  where  the  signal  is  30  dB  down  from  its  maximum 
value. 

If  one  allows  20  percent  of  the  channel  capacity  for  low  frequency  sensors 
(e.g.,  long-period  instruments,  etc.),  synchronization  and  framing,  error  coding, 
and  parity,  the  data  rate  is: 

R  =  120N  bits/second,  (F-l) 

where  N  is  the  number  of  short-period  sensors.  Selecting  the  2400  baud  line  as 
a  typical  voice  grade  line,  one  might  consider  the  smallest  instrument-converter- 
communications  component  group  as  20  short-period  instruments  sampled  at 
10  Hz  and  several  low  frequency  sensors  whose  aggregate  data  rate  is  less  than 
200  bits/second.  By  employing  many  identical  groups,  design  costs  are  mini¬ 
mized  and  a  direct  measure  of  sub-system  costs  becomes  available  in  terms  of 
a  system  performance  parameter. 

F.2.2  Gain 

Consider  an  array  of  sensors  in  which  it  is  desired  to  steer  many  conven- 

th 

tionally  *»2  formed  beams.  Let  x  (t)  be  the  output  of  the  n  sensor,  then  the 

th  ^ 

m  beam  output  is: 
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Figure  F-l.  Modified  Floating  Point— Data  Encoding  Modulo  4 
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where  N  is  the  number  of  sensors  in  the  array  and  Tmn  is  the  time  delay. 
Defining  gain  as  the  ratio  of  the  average  signal  power  to  the  average  noise 
power  of  the  beam  output  divided  by  the  ratio  of  the  average  signal  power  to  the 
average  noise  power  at  a  sensor,  it  can  be  shown  that: 


1+  (N-l)  P 
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where  pg  and  pn  are  the  average  signal  and  noise  correlation  coefficients, 
respectively. 

Equation  F-3  is  depicted  in  Figure  F-2  for  several  values  of  pn  as  a 
parameter.  First  note  that  only  when  p  is  very  small  does  the  denominator 


n 


of  Equation  F-3  approach  unity  (e.g.,  neglect  trivial  case).  On  the  other  hand, 


when  p  =  1/2,  a  loss  of  an  additional  3dB  is  incurred, 
s 


Using  the  classical  assumption,  that  is,  the  signal  is  correlated  (PR  =  1) 
and  the  noise  is  uncorrelated  (P  =  0),  Equation  F-3  reduces  to: 


G  =  10  l°g^Q  N  dB. 


(F-4) 


Unless  the  noise  correlation  is  small,  indiscriminant  addition  of  instruments 
will  not  always  yield  more  gain.  A  more  general  relationship  and  perhaps  one 
that  would  be  more  directly  applicable  to  array  synthesis  is  the  dependence  of 
the  correlation  coefficient  on  element  spacing. 

Let  us  further  assume  that  we  desire  to  place  a  large  number  of  elements 

in  a  given  space  such  that  we  are  interested  in  noise  correlation  coefficients  of 

_2 

the  order  of  10  .  Using  the  standard  deviation  as  an  indication  of  the  quality 

of  the  measurement,  it  can  be  shown  that: 
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where  W  is  the  noise  bandwidth  in  Hz  and  T  is  the  record  length  in  seconds. 
From  the  spectral  estimates  presented  in  Section  F.3.3,  we  select  W  =  1/2  Hz. 
If  cr^  is  to  be  less  than  the  value  of  p^,  we  require  a  minimum  record  length 
of  almost  three  hours. 

In  order  to  use  pair-by-pair  correlation  to  obtain  the  average  correlation 
coefficient,  the  effects  of  directional  noise  must  be  considered.  In  view  of  this 
and  the  aforementioned  restrictions,  it  seems  reasonable  to  suggest  that  corre¬ 
lation  measurements  be  compared  to  gain  measured  directly  by  conventional 
beamforming  and  average  power  techniques . 

F.2.3  Sampling 

In  Section  F.2.1,  it  is  suggested  that  a  sampling  rate  of  10  Hz  would  be 
adequate.  Since  data  reconstruction  techniques  permit  extension  of  the  effective 
sampling  rate,  high  frequency  data  can  be  retrieved  with  little  error.  (See 
Appendix  C.)  Since  this  procedure  would  not  be  applied  on  a  routine  basis  in  an 
operational  sense,  it  would  be  reasonable  to  sample  at  a  rate  commensurate 
with  the  real-time  requirement.  The  loss  in  conventional  processing  gain  as  a 
result  of  sampling1  is: 

5  log10e  2 

Gx  =  - (  2rf  ATy  dB,  (F-G) 

where  f  is  the  signal  frequency  in  Hz  and  AT  is  the  sample  period  in  seconds. 
Equation  F-6  is  depicted  in  Figure  F-3.  When  the  signal  frequency  is  less  than 
1.5  Hz,  a  loss  of  less  than  0.5 dB  is  incurred  when  the  sensors  are  sampled  at 
10  Hz. 

F.2.4  Time  Delays 

Conventional  beamformed  array  processing  requires  that  the  time  delays 
employed  in  Equation  F-2  be  determined  for  the  teleseismic  surveillance 
region  between  0.04  sec/km  and  0.08  sec/km. 
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Figure  F-3.  Sampling  Loss 
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These  delays  will  be,  in  general,  non-planar  functions  of  both  range  and 
azimuth,  and  it  is  planned  to  employ  existing  data  to  establish  the  initial  time 
delay  library.  Through  system  usage,  the  library  will  be  updated  and  improved. 
Interpolation  between  regions  for  which  data  exists  or  admittance  of  events 
whose  origin  is  not  in  the  teleseismic  region  may  reduce  the  number  of  calibra¬ 
tion  events  required.  Through  a  knowledge  of  event  hypocenter,  the  data  base 
for  the  inverse  velocity  space  to  geographic  coordinate  conversion  can  be 
developed. 

This  procedure  separates  array  phasing  and  event  location  into  two  tasks. 
Measurement  of  phase  speed  in  inverse  velocity  space  yields  the  appropriate 
time  delays  through  the  wave  front  approximations.  The  geographic  event  loca¬ 
tion  is  then  obtained  by  conversion  from  inverse  velocity  space  by  appropriate 
tables.  The  absolute  location  accuracy  is  dependent  upon  the  assumed  model, 
which  can  only  be  improved  through  calibration. 

An  error  in  the  time  delay,  therefore,  has  two  distinct  effects.  First,  the 
beam  is  defocused  and  a  loss3  is  incurred: 

G2  =  (10  log1Qe)  (27TfaT)2dB,  (F-7) 

where  f  is  the  signal  frequency  in  Hz  and  a-p  is  the  standard  deviation  of  the 
steering  error  in  seconds.  Equation  F-7  is  depicted  in  Figure  F-4. 

Referring  again  to  the  spectral  density  estimates  in  Section  F.3,  we  note 
that  little  signal  energy  is  received  above  1  Hz.  Using  the  1  Hz  curve  and 
limiting  the  loss  to  1  dB,  the  total  steering  (including  wave  front  approximations) 
standard  deviation  error  must  be  less  than  75  ms.  For  regions  where  this  is 
not  possible,  the  loss  manifests  itself  as  beam  defocusing  and  array  perform¬ 
ance  is  degraded.  It  is  also  pointed  out  that  the  steering  delays  error  may  be 
dependent  upon  array  aperture  and  no  attempt  has  been  made  here  to  take  this 
additional  loss  into  account. 

The  second  effect  is  the  uncertainty  in  event  location.  The  relative  ac¬ 
curacy  of  the  event  location  process  can  be  described  in  terms  of  the  system's 
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Figure  F-4.  Phasing  Loss 


ability  to  measure  time  in  a  noise  background  and  the  array  aperture.  Under 
Gaussian  assumptions  for  the  noise  background  and  statistical  independence 
of  noise  at  the  elements,  it  can  be  shown  that  the  root  mean  square  of  the 
relative  location  error  for  a  typical  range  angle  of  58°  (6400  Km)  is: 


7rfR„  2WT  S/N 

hi 


where  f  is  the  signal  frequency,  Rp>  is  the  radius  of  the  circular  array  and  WT 
is  the  gain-bandwidth  product. 

In  order  to  assess  the  absolute  location  accuracy  in  a  form  suitable  for 
experimental  evaluation,  Equation  F-8  should  be  extended  to  take  into  account 
signal  decorrelation  and  k-space  to  geographic  conversion  bias  effects.  It  is 
pointed  out  that  these  bias  errors  appear  experimentally  inseparable  and  are 
only  minimized  through  calibration.  Their  contribution  to  the  location  error 
has  not  been  taken  into  account  herein. 

F.  3  EXPERIMENTAL  RESULTS 

The  experimental  data  used  in  this  section  was  taken  from  LASA.  It  is 
presented  to  demonstrate  the  utilization  of  the  theory,  to  indicate  the  feasibility 
of  the  processing  concepts  hypothesized  in  Section  F.4,  and  to  depict  how  the 
results  of  field  tests  may  be  used  to  predict  expected  array-processing  system 
performance.  Although  the  data  sample  is  extremely  small,  it  is  believed  to 
be  typical  and  representative  of  prevailing  conditions  at  LASA. 

F.3.1  Gain 

For  the  particular  array  geometry  of  LASA,  several  of  the  subarrays  may 
be  combined  to  yield  many  different  array  configurations.  Using  selected  in¬ 
struments,  the  validity  of  the  classical  conventional  processing  gain  formula 
(10  log^N)  can  be  experimentally  verified  for  varying  conditions.  Depicted  in 
Figures  F-5  and  F-6  (solid  lines)  is  the  conventional  processing  loss  as  a 
function  of  the  number  of  LASA  sensors  utilized  for  three  array  apertures  (20, 
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Figure  F-5.  Conventional  Processing  Loss 
(0.7  -1.4  Hz) 
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Figure  F-6.  Conventional  Processing  Loss 
(0.9  -  1.4  Hz) 
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30,  and  50  kilometers).  In  Figure  F-5,  a  0.7-1. 4  Hz  third  order  Butterworth 
bandpass  filter  was  employed  and  in  Figure  F-6,  the  filter  3  dB  points  were 
0. 9-1.4  Hz.  The  dashed  lines  in  figures  F-5  and  F-6  (right  ordinate)  indicate 
the  number  of  sensors  divided  by  the  square  of  the  average  distance  between 
sensors  used  in  the  corresponding  gain  measurement.  The  scale  should  be 
modified  by  the  factor  l//r  when  comparing  the  quasi  density  expressed  here 
to  circular  area  densities. 

The  principal  conclusion  drawn  from  this  and  other  data  is  that  through 
proper  element  spacing,  conventional  processing1*2  and  appropriate  frequency 
filtering4,  it  is  possible  to  obtain  gains  which  are  within  about  1  or  2  dB  of 
theoretical  values  when  fifty  elements  are  densely  packed.  It  is  inferred  that 
through  proper  spacing,  considerably  more  elements  can  be  employed  with 
equivalent  results. 

F.3.2  Noise  Correlation 

The  average  element  spacing  available  from  LASA  required  to  support 
conventional  processing  is  perhaps  not  a  useful  design  parameter.  A  more 
conservative  criterion,  minimum  element  spacing,  will  also  be  considered. 

Exhibited  in  Figure  F-7  are  many  values  of  the  noise  correlation  coefficient 
as  a  function  of  the  spacing  between  instruments  for  one  sample  of  noise. 
Combinations  of  the  sensors  in  selected  LASA  subarrays  were  used;  and  the 
previously  specified  0. 9-1.4  Hz  filter  was  employed.  Note  that  the  average 
correlation  coefficient  passes  through  zero  at  approximately  five  kilometers. 
Because  of  the  anticipated  variations  in  the  zero  crossover  point,  we  will  employ 
a  modest  margin  and  use  6  kilometers.  Without  regard  to  beam  pattern,  a 
uniformly  distributed  array,  such  as  one  developed  by  contiguous  equilateral 
triangles  or  squares,  will  require  near  minimum  area  and  provide  near  maximum 
gain  for  a  small  number  of  sensors. 

In  order  to  employ  correlation  techniques,  the  noise  directionality  must  be 
taken  into  account.  To  support  the  inherent  isotropic  conjecture,  313  beams, 
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Figure  F-7.  Correlation  Coefficient  vs.  Distance  Between  Instruments 


whose  centers  formed  the  vertices  of  equilateral  triangles  of  leg  length  5  x  10  ** 

_2 

sec/Km  were  formed  in  the  circle  of  radius  8  x  10  sec/Km.  For  this  sample, 
the  maximum  deviation  of  each  beam  output  was  found  to  be  less  than  3  dB. 

The  temporal  stability  of  this  apparent  isotropic  condition  has  not  been 
investigated,  nor  has  a  statistically  significant  data  base  been  established  to 
indicate  that  these  results  are  typical. 

F.3.3  Spectral  Density 

The  pertinent  bulletin  characteristics  of  the  events4  presented  in  this 
section  are  tabulated  in  Table  F-l.  Their  spectral  density  functions  in  decibels 
referenced  to  one  nanometer,  are  exhibited  in  figures  F-8  and  F-9.  Typically, 
this  data  is  a  composite  of  13  seconds  of  event  data  and  50  seconds  of  record 
preceeding  the  event. 

For  Longshot  the  traces  peak  at  about  1.0  Hz.  Note  that  good  subarray 
and  array  gain  is  experienced  below  1.0  Hz.  However,  above  1.0  Hz,  the  sub¬ 
array  performs  almost  as  well  as  the  array.  Note  also  the  improvement  in 
subarray  gain  at  frequencies  above  1.0  Hz. 

TABLE  F-l 
BULLETIN  DATA 


Longshot 

Event 

Earthquake 

10/29/65 

Date 

4/8/66 

21:08:35.5 

LAO  Arrival 

01:56 :40 .2 

Rat  Island 

Location 

East  Coast  of  Kamchatka 

51.44°N  179.18°E 

Latitude/Longitude 

51.2°N  157 .7°E 

47.15° 

Range 

58.81° 

304.55° 

Azimuth 

313.0° 

6.1 

Magnitude 

5.9 
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Figure  F-8.  Longshot  S/N  Spectrum 
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Figure  F-9.  Kamctiatkci  S/N  Spectrum 


For  the  Kamchatka  event,  a  wider  signal  spectrum  is  observed;  however, 
the  principle  signal  spectrum  is  less  than  1.2  Hz.  Good  composite  subarray 
gain  is  achieved  throughout  the  spectrum,  while  good  subarray  gain  is  only 
achieved  above  1.0  Hz. 

Both  spectral  estimates  support  the  conjecture  that  noise  coherence  at  low 
frequencies  exists  and  from  Section  F.3.2  it  can  be  inferred  that  the  poor  sub¬ 
array  gain  is  due  to  the  element  spacing  within  the  subarrays. 

Since  the  signal  duration  of  Longshot  is  small  compared  to  the  record 
length  used,  it  is  suggested  that  the  poor  composite  subarray  performance 
above  1  Hz  may  be  due  to  signal  reverberation.  This  has  not  been  experimentally 
tested,  nor  has  there  been  sufficient  data  analyzed  to  warrant  this  implication. 


F.4  SYSTEM  HYPOTHESIS 

The  array-processing  requirements  naturally  divide  into  detection  and 
event  functions.  1>2  The  detection  function3  primarily  determines  the  existence 
of  an  arrival  and  its  approximate  location.  The  association  of  several  arrivals 
into  an  event  and  the  estimation  of  event  magnitude,  location  and  depth  is  the 
event  function.4  Classification  is  considered  a  post-event  processing  function 
and  beyond  the  scope  of  this  document. 

The  detection  function  is  primarily  a  real-time  function  (i.e.,  it  must  be 
accomplished  in  a  manner  so  as  to  not  fall  behind),  while  the  event  function 
primarily  requires  a  processing  capability  based  on  the  event  occurrence  rate. 
The  detection  function  operates  on  a  power  threshold  ratio  and  decision  logic 
criterion  over  a  frequency  band  where  the  signal  to  noise  ratio  peaks.  The 
event  function,  on  the  other  hand,  requires  maximal  gain  and  fidelity  (or  at 
least  more  than  the  detection  function)  to  insure  waveform  preservation, 
arrival  identification,  location  accuracy  and  an  appropriate  evasion  capability 
level.  The  detection  function  implies  area  surveillance.  To  insure  arrival 
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inclusion  in  event  processing,  the  event  function  requires  only  a  beam  coverage 
density  over  the  area  of  the  detection  location  uncertainty. 

In  addition,  the  detection  function  is  normally  specified  by  a  detection  capa¬ 
bility  and  a  false  alarm  probability,  while  the  event  function  is  specified  by 
geophysical  parametric  estimates  of  arrival  magnitude,  hypocenter,  and  time  of 
origin.  Typically,  more  information  is  required  in  event  processing  than  for 
detection  processing. 

The  above  noted  dichotomies  strongly  suggest  a  multi-level  array  con¬ 
figuration.  Depicted  in  Figure  F-10  is  a  hypothetical  array-processing  system 
from  which  some  quantitative  performance  trade-offs  can  be  estimated. 

To  provide  a  basis  for  design  synthesis,  the  system  objectives  require 
specification.  For  detection,  we  will  arbitrarily  require  a  conventionally 
beamformed  array  gain  of  20dB.  The  detection  capability  can  be  expressed  as: 

Dc=VTc’dB’  <F-9> 

where  Gp  is  the  total  processing  gain,  including  both  filter  and  array  gain  and  T 
is  the  detector  threshold  margin.  For  reasonable  false  alarm  rates  anticipated 
from  the  threshold  decision  logic4,  we  conservatively  assume  that  the  filter  gain 
and  the  threshold  margin,  are  approximately  equal.  Then,  Dc  is  essentially  the 
gain  with  respect  to  a  filtered  trace;  and,  therefore,  N=100. 

For  the  event  function  objectives  considered  here,  the  location  requirement 
is  probably  the  most  demanding.  The  event  array  capability  can  be  expressed 
in  terms  of  the  event  array  gain  and  its  filter  gain,  namely: 

EC  =GE  +GF’  dB>  (F_10) 

However,  the  event  location  is  a  function  of  both  event  array  gain  (neglecting 
filter  gain),  and  the  array  aperture  expressed  as  the  event  array  radius  (R^). 
Here,  we  will  arbitrarily  require  a  location  accuracy  of  50  kilometers  when  the 
signal  to  noise  ratio  (S/N)  on  the  processed  output  is  about  12  dB. 
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Figure  F-10.  System  Parameters 
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In  the  following  paragraphs,  the  system  performance  requirements  will  be 
developed  in  terms  of  the  array  specification  parameters. 

F .  4 . 1  Subarray  C  ommunications 

Using  the  grouping  of  sensors  suggested  in  Section  F.2.1  and  F.3.2,  con¬ 
sider  an  equilateral  triangular  arrangement  of  instruments  as  shown  in  Figure 
F-ll.  Denoting  the  minimum  distance  between  elements  as  r,  the  area  of  the 
circumscribed  circle  containing  N  instruments  is  represented  by  its  radius. 
Using  r  =  6  kilometers  [to  obtain  zero  noise  correlation  (Section  F.3.2)],  we  note 
from  Figure  F-12  that  the  minimum  grouping  radius  (R)  for  twenty  sensors  is 
about  12.5  kilometers.  Although  crude,  we  are  able  to  suggest  that  nineteen 
sensors,  when  placed  on  concentric  hexagons  (and  one  in  the  center)  with  a 
minimum  element  separation  of  6  kilometers,  will  occupy  an  area  of  radius 
12.5  kilometers,  yield  a  gain  of  almost  13  dB  when  processed  with  conventional 
beamforming  techniques  and  can  be  supported  by  a  2400  baud  line.  In  the 
following,  a  subarray  aperture  (2R^)  of  25  kilometers  wall  be  employed. 

F.4.2  Detection  Array-Processing 

Assuming  a  uniform  equilateral  distribution  of  sensors,  the  number  of 
beams  required  to  cover  the  teleseismic  region  between  0.04  sec /kilometer 
and  0.08  sec/kilometer,  with  a  loss  of  no  more  than  3dB,  can  be  related  to 
detection  array  radius  (Rq)  and  the  frequency  of  the  received  signal,  as  ex¬ 
hibited  in  Figure  F-13. 

Using  Figure  F-12,  we  require  the  array  radius  to  be  about  30  kilometers 

to  achieve  20dB  array  gain,  when  the  minimum  spacing  between  elements  is 

6  kilometers.  For  a  1  Hz  signal  and  an  array  aperture  of  60  kilometers,  we 

t 

note  from  Figure  F-13,  that  about  100  beams  are  required  to  cover  the 
teleseismic  region. 

The  detection  processing  (DP)  is  thereby  grossly  specified.  To  wit,  apply 
one-step  beamforming  (BF),  filtering  (FIL)  and  rectification-integration  and 
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Figure  F-ll.  Equilateral  Triangle  Distribution 


F-22 


(9P)  NIV9 


F-23 


Figure  F-12.  Detection  Array  Gain 
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Figure  F-13.  Teleseismic  World  Coverage 
(0.  04-0.  08  Sec/Km) 
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thresholding  (RIT).  Beamforming  and  filtering  (3rd  order  bandpass)  will  be 
executed  at  10  Hz  and  RIT  at  2.5  Hz.  The  output  will  be  a  queue  of  detected 
arrivals  containing  identification,  the  time  which  the  detection  criterion  was 
first  satisfied  and  first  not  satisfied,  maximum  magnitude  per  beam  satisfying 
the  criterion  and  the  average  noise  when  the  arrival  was  declared. 

To  assess  the  requirement  relative  to  existing  capability4,  a  comparison 
to  the  micro-programmed  machines  was  made,  and  is  depicted  in  Figure  F-14. 
Using  100  sensors  uniformly  placed  over  an  aperture  of  60  kilometers,  the 
speed  and  memory  requirements  are  plotted  as  a  function  of  the  number  of 
beams  formed.  When  generating  100  detection  beams  from  100  sensors,  the 
memory  requirements  are  well  within  that  available;  and  the  speed  requirement 
is  about  70  percent  of  present  capability. 

F.4.3  Event  Array-Processing 

Application  of  the  theory  presented  in  Section  F.2.3  implies  that  the  signal 
to  noise  ratio  of  the  fully  processed  event  record  depends  upon  event  magnitude. 
For  example,  a  magnitude  4.0  event  will  be  some  10 dB  stronger  than  a  magni¬ 
tude  3.5  event  and  one  may  expect  the  output  record  signal  to  noise  ratio  (S/N) 
to  reflect  this  difference.  It  is,  therefore,  reasonable  to  specify  the  location 
precision  in  terms  of  not  only  the  event  array  aperture,  signal  frequency  and 
system  gain-bandwidth  product,  but  the  event  processing  output  S/N  as  well. 

Equation  F-8  is  depicted  in  Figure  F-15.  For  this  example,  we  have  taken 
f  =  1  Hz  and  WT  =  1.  If  one  desires  an  event  location  error  of  about  50  kilo¬ 
meters  when  S/N  =  12  dB,  a  200  kilometer  event  array  aperture  is  required. 

In  order  to  support  this  requirement,  the  event  beams  are  specified  in 
terms  of  the  uncertainty  of  event  location  as  determined  by  the  detection 
processing  function.  In  Figure  F-16,  the  number  of  event  beams  is  presented 
as  a  function  of  the  ratio  of  the  event  and  detection  array  radii  with  received 
signal  frequency  as  a  parameter.  The  9dB  detection  beam  contour  is  a  con¬ 
servative  selection.  To  demonstrate  sensitivity,  the  3  dB  contour  is  also 
depicted.  Referring  to  Section  F. 3.3,  there  is  very  little  energy  beyond  2  Hz. 
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Figure  F-14.  Detection  Processing 
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At  the  3  dB-3  Hz  point,  we  take  Rg/Rj)  =  3  and  no  more  than  200  event  beams  are 
required. 

The  remaining  parameter  needed  to  specify  the  event  beamforming  require¬ 
ment  is  the  number  of  sensors  (or  groups)  constituting  the  event  array.  For 
convenience,  we  will  assume  that  the  number  of  event  sensors  (Nj?)  is  equal  to 
the  number  of  detection  sensors  (Nd).  The  effect  of  more  or  less  sensors  will 
be  examined  subsequently. 

Defining  five  subarrays  of  19  or  20  sensors  each,  execute  two-step  beam¬ 
forming  generating  event  array  beams  from  one  subarray  beam  per  subarray. 
The  subarray  beams  and  filters  are  assumed  to  be  formed  at  a  rate  of  10  Hz 
with  a  beamforming  resolution  of  10  Hz.  The  array  beam  resolution  is  10  Hz, 
but  is  formed  at  5  Hz;  and  RIT  is  executed  at  2.5  Hz.  The  machine  speed  and 
memory  requirements  are  depicted  in  Figure  F-17.  At  200  beams,  the  speed 
is  less  than  20  percent  of  available  real-time  capacity  and  the  memory  is  much 
less  than  the  detection  processor. 

It  would,  however,  be  prudent  to  use  both  the  detection  array  and  the  event 
array  instruments  for  event  processing  (EP).  In  this  manner,  an  additional 
theoretical  gain  of  3  dB  is  potentially  available.  Here,  the  processing  require¬ 
ment  is  about  equivalent  to  available  real-time  capability,  and  the  memory  re¬ 
quirement  is  essentially  the  sum  of  the  detection  and  event  processing  needs. 
Since  event  processing  speed  capability  is  related  to  the  rate  of  event  occur¬ 
rence,  overall  event  beam  forming  computations  are  deemed  small  compared 
to  the  arrival  and  best  beam  selection  process,  and  the  bulletin  computations. 
These  are  the  subject  of  current  study. 

F.4.4  Event  Location  Confidence 

The  relative  location  accuracy  prediction,  presented  in  Section  F.4.3,  is  a 
bias  free  statistic  and  subject  to  a  probability  confidence.  The  confidence  error 
sensitivity  bound  is  depicted  in  Figure  F-18.  If  location  estimates  are  to  be 
made  with  greater  than  50  percent  confidence,  select  a  =  1.0. 
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Figure  F-l  7.  Event  Processing 
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Figure  F-18.  Minimax  Sensitivity  Bound 
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To  relate  the  process  to  its  input,  the  number  of  detection  and  event  sensors 
has  also  been  introduced.  In  Equation  F-8,  the  product  of  the  circular  event 
array  radius  and  the  ratio  of  the  root  of  the  number  of  event  and  detection 
sensors  has  been  modified  by  the  confidence  error  sensitivity  parameter  (i.e. , 

Rg-yl+Ng/N^1  /  100a).  Under  the  assumption  there  is  no  signal  loss  through 

the  process,  system  performance  can  be  estimated  in  terms  of  event  magnitude 
with  respect  to  the  input  ambient  noise  conditions  as  depicted  in  Figure  F-19. 

To  estimate  the  number  of  event  sensors  required  to  achieve  the  event 
location  objective,  specification  of  the  input  event  magnitude  is  necessary.  For 
example,  we  will  require  with  50  percent  confidence,  a  100  kilometer  event 
location  accuracy  for  events  up  to  one  magnitude  unit  below  the  ambient  noise 
level.  For  the  previously  discussed  case  (Ng  =  Np),  we  note  from  Figure  F-19 
that  Rj?  rs160  kilometers.  However,  if  Np/Npj  =  3,  the  event  array  radius  can 
only  be  reduced  to  about  110  kilometers. 

It  is  reiterated  that  both  the  steering  delay  and  geographic  coordinate  con¬ 
version  errors  have  been  neglected.  Inclusion  of  the  effect  of  these  bias  errors 
in  this  random  variable  portrayal  will  no  doubt  worsen  the  location  error 
estimates. 

F.4.5  Beam  Patterns 

The  preceding  discussion  has  been  without  regard  to  array  geometry  effect 
on  beam  patterns.  Depicted  in  Figure  F-20  is  the  beam  pattern  of  a  21  element 
quasi  uniformly  distributed  pentagon  25  kilometer  aperture  subarray.  No  side 
lobes  greater  than  about  9.0  dB  are  noted  out  to  0.2  sec /kilometer.  A  60  kilo¬ 
meter  array  with  a  similar  distribution  and  with  about  100  elements  can  be 
synthesized  and  when  its  subarray  beam  is  combined  with  the  subarray  beams 
of  the  25  kilometer  aperture  subarrays,  the  composite  array  side  lobe  per¬ 
formance  will  be  somewhat  improved. 

Depicted  in  Figure  F-21  is  the  beam  pattern  of  a  regular  pentagon  event 
array  with  one  element  placed  in  the  center  for  a  total  of  six  elements  (subarrays). 
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Figure  F-19.  Event  Array  Parameter 
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Figure  F-20,  Quasi  Pentagon  Subarray  Beam  Pattern 
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Figure  F- 21 .  Regular  Pentagon  Beam  Pattern 


The  radius  is  taken  as  100  kilometers  which  results  in  a  leg  length  of  about 
85  kilometers.  The  first  side  lobes  exist  at  0.01  sec/kilometer  and  are  only  3  dB 
down.  The  dashed  concentric  circles  labeled  (1)  and  (2)  are  the  subarray  beam 
1  and  2  dB  contours  respectively.  Due  to  the  subarray  beams,  the  event  array 
side  lobe  will  be  down  at  least  4.5  dB  at  0.01  sec/kilometer.  The  technique 
employed  here,  using  the  subarray  rejection  capability  to  reduce  event  array 
side  lobes,  can  be  augmented  by  additional  element  placement. 

When  specifying  the  lumped  detection  array  and  a  distributed  event  array 
primarily  for  the  location  requirement  (and  for  the  long-period  instrument 
spacing),  it  is  desirable  to  obtain  maximum  utilization  of  the  sensors  in  both 
arrays.  As  suggested  in  Section  F.4,  the  detection  and  event  array  require¬ 
ments  are  vastly  different  and  support  different  specifications.  Hence,  a 
compromise  is  needed. 

Consider  the  case  when  Ng  =  Np>,  and  we  place  the  detection  array  in  the 
center  of  an  event  array.  In  this  case,  one-half  of  the  total  sensors  do  not 
materially  contribute  to  event  location.  However,  when  the  detection  array  is 
placed  on  peripheral  of  the  event  array  (as  shown  in  Figure  F-10),  most  of  the 
elements  contribute  to  event  location.  This  argument  supports  an  asymmetric 
array  geometry  and  identifies  the  as  yet  to  be  investigated  array  azimuthally 
orientated  noise  directionality  effects. 

The  asymmetrical  array  loss,  estimated  by  element  weighting,  was  found 
to  be  less  than  0.4dB  for  values  of  Ng/Nj)  =  2,  and  negligible  for  Ng/Nj)  <  2. 

The  asymmetrical  array  beam  pattern  shown  in  Figure  F-22  was  computed  by 
weighting  the  six  element  case  above  such  that  Np>  =  Ng.  Here,  there  appears 
significant  degradation  in  side  lobe  structure  and  further  array  design  study  is 
required  before  an  event  array  can  be  hypothesized.  However,  the  effects  of 
array  aperture  on  identification  of  post-P  arrivals  and  evasion,  and  the 
presently  envisioned  classification  procedures  support  the  above  suggested 
gross  parametric  relations.  Array  shading  can  be  introduced  and/or  additional 
subarrays  may  be  placed  contingent  upon  long  period  instrument  requirements 
to  minimize  side  lobes  and  directional  noise,  and  maximize  location  accuracy 
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Figure  F-22.  Asymmetrical  Pentagon  Beam  Pattern 


along  certain  azimuths.  The  potency  and  practicality  of  these  suggestions  have 
not  been  investigated,  nor  has  the  effect  of  time  delay  errors  as  a  function  of 
array  radius  been  established. 

F .  5  SUMMARY 

The  tri-level  array  processing  system  hypothesized  has  been  specified  in 
terms  of  gross  system  requirements.  The  pertinent  system  parameters  are 
summarized  in  Table  F-2. 

This  analysis  has  been  restricted  to  consider  only  a  portion  of  the  array¬ 
processing  design  problem  in  terms  of  a  particular  desired  product.  In  order 
to  demonstrate  the  approach,  certain  objectives  were  specified  (i.e.,  20  dB 
detection  array  gain  and  50  kilometer  relative  location  accuracy  when  the  output 
record  S/N  is  greater  than  12  dB).  Other  products  may  be  desirable.  For 
example,  integration  of  the  long-period  instruments,  although  provided  for  in 
the  subarray  communication  loading  and  event  array  geometry,  is  a  glaring 
omission.  It  can,  however,  be  readily  incorporated  with  a  knowledge  of  the 
long-period  noise  coherence.  Further,  an  isotropic  noise  field  has  been  assumed, 
and,  for  site  locations  where  a  strong  directional  noise  source  is  likely  to  exist, 
this  conventional  processing  synthesis  procedure  may  require  extension. 

Further,  the  effects  of  detection  threshold  margin  on  false  alarm  rate,  the 
effects  of  aperture  on  steering  delay  errors,  and  the  effects  of  the  apparently 
inseparable  steering  delay  and  geographical  bias  errors  on  location  accuracy 
have,  in  the  interest  of  expediency,  been  ignored.  When  suitably  accurate  ex¬ 
perimental  data  recorded  and  processed  with  sufficient  precision,  becomes 
available,  it  is  planned  to  extend  this  analysis  as  needed. 

Using  the  grouping  concept  hypothesized  in  Section  F.4.1,  for  both  the  de¬ 
tection  and  event  arrays,  we  note  that  we  are  able  to  withstand  subarray  failures 
without  catastrophic  system  effects.  If  the  redundant  processor  concept2*3  is 
maintained  for  detection  processing  and  sensor  recording,  the  event  processor 
may  assume  the  detection  function  for  the  case  of  a  detection  channel  failure. 

The  additional  requirement  imposed  is  one  of  sizing  the  event  processor  such 
that  it  has  a  sufficient  catch-up  capability. 
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Table  F-2 


PARAMETER  SUMMARY 


Case 

Objective 

Parameter 

Subarray 

2400  Baud  Line 

Uniform  Distribution:  20  SP 

1  LP 

78  dB  Range 

Subarray  Aperture:  25  Km 

1/4  dB  Resolution 

Sampling:  10  Hz  (5  Hz  Filter)  SP 

1  Hz  (0.5  Hz  Filter  )  LP 

Detection  Array 

20  dB  Detection  Gain 

Uniform  Distribution:  100  SP  or  5  SA 

0.04  -  0.  08  sec/Km 

Array  Aperture:  60  Km  or  5  SA 

3  dB  Beam  at  1 .5  Hz 

Beams:  100 

Processing:  Conventional  (One -Step) 

Event  Array 

50  Km  Relative  Location 

Circular  Distribution:  5  Subarrays 

0.5  Mag.  Event  Below  Noise 

Array  Aperture:  200  Km 

Confidence  of  50% 

Beams:  1  Subarray  Beam  and 

200  Array  Beams 

12  dB  Processed  Record  S/N 

Processing:  Conventional  (Two-Step) 
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With  respect  to  minimizing  installation  and  operating  cost,  the  above 
delineated  synthesis  procedure  provides  a  cost  function  which  is  strongly 
dependent  on  array  radius  and  the  total  number  of  sensors.  The  constraint  is 
related  to  both  system  objectives. 
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